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DIGITAL-COMPUTER PROGRAM FOR DESIGN -ANALYSIS OF 
SALIENT, WOUND POLE ALTERNATORS 
by David S. Repas 
Lewis Research Center 

SUMMARY 

A digital computer program for analyzing the electromagnetic design of salient, 
wound pole alternators, is presented. The program, which is written in the FORTRAN IV 
programming language, calculates the open-circuit saturation curve, the field-current 
requirements at rated voltage for various loads and losses, efficiency, reactances, 
time constants, and weights. The methods used to calculate some of these items are 
presented or appropriate references are cited. 

The program is useful for parametric studies and alternator design optimization. 
Instructions for using the program and typical program input and output for an alternator 
design are given in the appendixes. Also included are an alphabetical list of FORTRAN 
symbols and the complete program listing with flow charts. 


INTRODUCTION 

Dynamic energy conversion systems developed at the NASA Lewis Research Center 
for use in space have used solid rotor, brushless alternators. As a result, computer 
programs for analyzing the homopolar inductor and the outside-coil Lundell alternators 
were developed. These programs are described in references 1 and 2. 

In recent months, researchers at Lewis have been studying the concept of an inte- 
grated engine generator for aircraft secondary power systems (ref. 3). For this appli- 
cation a salient, wound pole alternator was considered. The salient, wound pole alter- 
nator is inherently much lighter in weight than a solid-rotor alternator and usually has 
lower reactances resulting in better voltage regulation and transient performance. 

Because of these considerations, it became necessary to develop a salient, wound 
pole computer program similar to the solid-rotor alternator programs. The fastest and 
easiest way to develop such a program was to modify the existing Lundell program. In 



addition, the material presented in reference 4 was used extensively, particularly in 
calculating reactances and time constants. 

Most of the necessary modifications arise directly from the differences in the mag- 
netic circuits of the salient, wound pole and Lund ell alternators. Also the effect of rotor 
damper bars on alternator performance was included in the salient, wound pole program. 
While the resultant salient, wound pole alternator program differs from the Lundell 
computer program, there still remain many similarities. Both programs are written 
in FORTRAN IV programming language, most of the FORTRAN symbols used are the 
same, and the input requirements and output formats are similar. 

This report describes the program in detail and gives some of the methods of calcu- 
lation. Instructions for using the program, typical input and output, the program listing 
and flow charts, and a FORTRAN symbol list are included in the appendixes. 


COMPUTER PROGRAM DESCRIPTION 
General Description 

The salient, wound pole alternator computer program is an analysis program. This 
means that the program accepts as input a complete electromagnetic alternator design; 
from this, it calculates the open-circuit saturation curve, the field-current requirements 
at rated voltage for various loads and losses, efficiency, several reactances and time 
constants, and weights of electromagnetic components. The results of the calculations, 
together with the input, are then printed out to provide a complete, self-explanatory rec- 
ord. An explanation of how to use the program, including a list of input variables with 
definitions, is given in appendix A. A typical program output is shown in appendix B. 

The program is written in FORTRAN IV and requires approximately 14 200 storage 
locations for execution. At Lewis, the program has been used on the IBM 7044-7094 
Mod II direct-couple system using a FORTRAN IV version 13 compiler. For this system 
typical pre-execution time is 1. 10 minutes and typical execution time is 0. 04 minute per 
alternator design. 

Figure 1 shows a simplified flow chart of the complete program. As can be seen, 
the program consists of a main program SALENT and four subrountines SUBSAL, 
OUTPUT, MAGNET, and WIND GE. The subroutines are necessary, in part, because 
one program is too large to compile with the available core storage locations. FOR- 
TRAN listings and detailed flow charts are given in appendix C. 

Communication between the main program and its subroutines is by means of COM- 
MON except for subroutine WINDGE which uses an argument list. 
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Description of Alternator to Which Program is Applicable 


The basic alternator for which the program was written is the salient, wound pole 
alternator. A cross-sectional view of this type of alternator showing the main flux path 
is given in figure 2. For clarity, a four-pole alternator is shown. In addition to the 
main flux, leakage fluxes between the rotor poles are present. 

The main flux flows from a rotor north pole, across the air gap and then radially 
through the stator teeth into the stator back-iron or core. It then goes circumferentially 
through the stator core for a distance of one pole pitch, enters the teeth and continues 
radially across the air gap into a rotor south pole. From there the flux goes radially 
through the south pole, circumferentially through the rotor back-iron or core and then 
completes its path radially through a north pole. 

A number of assumptions, in addition to those implicit in the geometric configura- 
tion, are made regarding the alternator: 

(1) Rotor poles and core are made of the same material. 

(2) Alternator armature winding is three-phase and Y connected. 

(3) The field coil cross-section is rectangular in shape. 

(4) There is one field coil per pole. 

(5) The eddy current loss in round conductors is assumed to be zero. 

(6) The only leakage flux paths are between the rotor north and south poles. 

In contrast to the restrictions imposed on the alternator by the preceding assumptions, 
several options are available to the program user. These options, which increase the 
applicability of the program, are as follows: 

(1) Armature conductors maybe round or rectangular. 

(2) Field conductors may be round or rectangular. 

(3) Armature conductors may be stranded. 

(4) Rotor and stator need not be made of the same magnetic material. 

(5) Damper windings may or may not be present. If there is no damper winding, the 
rotor pole face is assumed to be smooth, that is, without slots. 

(6) If the damper windings are present, the damper bars may be either round or 
rectangular. 

(7) Five different stator slot configurations may be used. 


METHOD OF CALCULATION 

This section of the report outlines in general terms the method of calculation used 
in the computer program. However, because of the length of the program and the large 
number of equations involved, specific equations will not, except in a few instances, be 
given. These equations can be found in volume 2 of reference 4. In addition, detailed 
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information and equations may be found in the program listings and flow charts in ap- 
pendix C. To assist in locating information in the listing COMMENT, cards are used to 
identify the major calculations. Of further value is appendix D which is a FORTRAN 
symbol list. 

Many variables, both in the report and in the computer program, are expressed in 
the per unit system. Per unit quantities are defined as follows: 

1 per unit voltage = rated voltage 

1 per unit current = rated current 

1 per unit volt-amperes = rated volt-amperes 

1 per unit impedance = (line-to-neutral rated voltage)/(rated current) 

Some variables are given in percent; these are per unit values multiplied by 100. Equa- 
tions given throughout this report will use the same FORTRAN symbols used in the pro- 
gram and as defined in appendix D. Where a FORTRAN symbol does not exist, an or- 
dinary algebraic symbol will be used. 


Leakage Permeance Calculations 

The rotor leakage flux has a pronounced effect on the saturation curves and, to a 
lesser degree, on alternator efficiency and weight. Because the leakage paths (fig. 3) 
are through air (or some other gas exhibiting linear magnetic properties), they can be 
characterized by a permeance value. The permeance of a leakage path is defined as 
the ratio of flux through the path to the ampere-turn drop across the path. 

As indicated in figure 3, four rotor leakage permeances are considered in the pro- 
gram. These are 

(1) PI - permeance between the inner pole head surfaces 

(2) P2 - permeance between the pole head end surfaces 

(3) P3 - permeance between the inner pole body surfaces 

(4) P4 - permeance between the pole body end surfaces. 

These permeances are calculated using the equations presented in reference 5 (pp. 73 
to 75 and 273 and 274). 


Magnetic Calculations 

The purpose of the magnetic calculations is to compute the flux densities throughout 
the alternator, the ampere-turn drop across the various parts of the magnetic circuit, 
and the total field excitation. In the program the computations are divided into three 

categories: 
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(1) No-load rated-voltage magnetization characteristics 

(2) Alternator load characteristics - This category comprises magnetic calculations 
at rated voltage for several loads at a fixed power factor. These calculations may then 
be repeated for up to five different power factors. 

(3) No-load saturation data - This category gives up to 10 points on the open-circuit 
saturation curve starting at voltage VMIN. 

The magnetic calculations are carried out, for the most part, in subroutine MAG- 
NET, which is called by the main program SALENT as required (see fig. 1). However, 
before calling MAGNET for the first time, all areas and lengths needed by MAGNET are 
calculated in SALENT. 

Just before each transfer to the subroutine, three additional quantities needed by 
MAGNET are calculated in SALENT. These quantities are the useful flux per pole in the 
air gap (PPL), the ampere-turn drop in the air gap (FGL1), and the demagnetizing 
ampere-turns due to armature current (FGXL). The quantities PPL, FGL1, and FGXL 
are functions of voltage, load, and load power factor. 

At no load, PPL and FGL1 are proportional to the output voltage and, since the 
armature current is zero, the demagnetizing ampere-turns FGXL are zero. Expressed 
in equation form: 


PPL = (EDD)(FQ) 

FGL1 = (EDD)(FH) 

FGXL = 0 

where FQ is the useful flux per pole at no-load rated voltage (calculated in SALENT), 
FH is the air-gap ampere-turn drop per pole at no-load rated voltage (calculated in 
SALENT), and EDD is the constant of proportionality, here, equal to per unit alterna- 
tor output voltage. 

Under load (at rated voltage), the calculation of PPL, FGLl, and FGXL requires 
reference to the alternator phasor diagram shown in figure 4 (ref. 6, p. 238). Con- 
struction of the diagram starts with the terminal voltage and armature current phasors 
which, along with the angle between them, are determined by the particular load for 
which PPL, FGLl, and FGXL are to be calculated. The phasors representing resis- 
tance and reactance drops are then added to the diagram using values of armature resis- 
tance and synchronous reactance values previously calculated in SALENT. This locates 
the quadrature axis. 

The angle between the quadrature axis and the current phasor is of interest for cal- 
culating FGXL. This angle, labeled XI and called the internal power factor angle, can 
be calculated from the geometry of the phasor diagram. The length OP is also deter - 
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mined from the phasor diagram. The length OP is the quadrature axis air-gap voltage 
modified to account for flux wave form distortion under load (fig. 4). The quantities 
PPL, FGL1, and FGXL can now be determined from the following equations: 

FGXL = FGML • G • sin(Xl) 


For power factors <0. 95 


PPL = FQ • (length of OP) 
FGLl = FH • (length of OP) 


for power factors >0. 95 


PPL = FQ • (length of OP) • 1. 10 
FGLl - FH • (length of OP) • 1. 10 

where FGML is the demagnetizing ampere-turns at rated load and G is the volt- 
ampere output of the alternator (in per unit). The PPL, FGLl, and FGXL calculations 
are summarized in table I. 

With PPL, FGLl, and FGXL known, the calculations are turned over to subroutine 
MAGNET. The approach taken in MAGNET is to represent the actual magnetic circuit 
of the alternator by an equivalent electrical circuit having lumped resistive elements. 

The resistors representing air gaps or leakage flux paths are assumed linear. Those 
that represent a part of the flux path through iron must be considered nonlinear. The 
equivalent circuit over one pole pitch for the salient, wound pole alternator is given in 
figure 5. 

From the useful flux per pole (PPL) the flux densities and ampere-turn drops for the 
stator teeth (FTL) and stator core (FCL) are determined. These ampere-turn drops, 
together with the ampere-turn drop in the air-gap (FGLl) and the demagnetizing ampere- 
turns (FGXL), are used to determine the rotor leakage flux (PPRL). With the magneto- 
motive force (mmf) acting on each leakage flux path known, the leakage flux can be cal- 
culated from 


Leakage flux = mmf • permeance 

The mmf across paths 1 and 2 (fig. 3) is 

2 • (FGLl + FGXL + FTL + FCL) = 2 • mmf tot 
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For paths 3 and 4, the mmf is 2 • mmfj. t at the pole head and zero at the pole base, the 
average value being equal to mmf^. The rotor leakage flux is added to the useful flux 
to obtain total flux. The total flux is used to compute the flux densities and ampere-turn 
drops for the rotor pole and rotor core. This completes the equivalent circuit analysis. 


Efficiency and Loss Calculations 

Individual losses and efficiency are calculated at rated voltage for several loads of 
increasing magnitude, continuing until the alternator saturates or until calculations have 
been completed for five loads. While the first load at which loss calculations are made 
must always be zero, the program user has the option of specifying any or all of the re- 
maining four loads. These loads are designated by G within the program (G is the 
volt-ampere alternator output expressed in per unit or percent). 

These calculations are first carried out for rated power factor. They are then re- 
peated up to four times for any other power factor specified in the program input. The 
individual losses that are calculated by the program, along with the method of calcula- 
tion or references, are presented in the following sections. 

Field conductor losses (PR) and armature conductor losses (PS) . - These losses are 
given by the expression I 2 R where I is the direct- current (dc) or root-mean- square 
(rms) current in the winding, as appropriate, and R is the direct current winding resis- 
tance corrected for the winding temperature. Correcting the winding resistance for 
temperature involves several assumptions: 

(1) The average no-load temperature is known or can be estimated. 

(2) The average rated-load winding temperature is known or can be estimated. 

(3) The average winding temperature is a parabolic function of the current in the 
winding. 

The assumed variation in the armature winding temperature as a function of current 
is illustrated in figure 6(a). A similar curve for the field winding is shown in figure 6(b). 

Armature conductor eddy loss (EX) . - References 4 and 7 present a discussion of 
armature conductor eddy losses. In the program these losses are assumed to be zero 
for round conductors because round conductors are usually small. For rectangular con- 
ductors they are given by 


EX = (EZ - 1)(PS) — 

HM 

where EZ is the eddy factor, which is a function of conductor geometry, armature 
winding parameters, frequency, and armature conductor resistivity; PS is the armature 
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conductor loss; CL is the length of armature conductor within the stator slot; and HM 
is the total armature conductor length. 

Pole face losses (WN and PP) . - The no-load pole-face loss (WN) (see refs. 4 
and 8) is calculated from the equation 

WN = (GA)(D1)(D2)(D3)(D4)(D5)(D6) 

where GA is the main air gap area and D1 to D6 are empirical factors that are func- 
tions of, respectively, pole-face lamination thickness, air-gap flux density, slot fre- 
quency, slot pitch, slot opening to air-gap ratio, and pole embrace. As given in refer- 
ence 4 the empirical curve that defines the factor D3 is limited to a maximum slot 
frequency of 5800 hertz. In the program, for all slot frequencies greater than 5800 
hertz, values of D3 are obtained by extrapolation. 

Under load (ref. 9, eq. (22)) the pole face loss (PP) is given by 

PP = (GM)(WN) GM > 1 

where GM is the function of alternator current, ratio of slot opening to air gap, number 
of conductors per slot, number of parallel circuits, and air-gap ampere-turns at no-load 
and rated voltage. 

Damper losses (WW and DL) . - The no-load damper bar loss (WW) is calculated as 
shown in reference 10 using the cold damper-bar temperature (T33). Under load (ref. 9, 
eq. (22)) the damper bar loss (DL) is given by 

DL = (GM)(WW) GM > 1 

For this calculation, the hot damper bar temperature (T3) is used regardless of the 
magnitude of the load. 

Stator core loss (WQL and WQ) and stator tooth loss (ST and WT) . - The respec- 
tive equations used to calculate these losses are 


stator core loss = 3(stator core weight)(WL) 


stator tooth flux density V 
BK / 


stator tooth loss = 3(stator tooth weight)( WL) ( ^ j ^Jo ^ h nux density \ 

V BK / 

where WL is the core loss at flux density BK and at rated alternator frequency and 
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BK is the flux density at which WL is measured. Both WL and BK are obtained 
from the data provided by the material manufacturer. 

Miscellaneous load loss (WMIS) . - In an alternator there are generally additional 
electromagnetic losses not accounted for by the losses enumerated so far. These losses, 
labeled miscellaneous load loss, are assumed to be 1 percent of the kilovolt-ampere out- 
put of the alternator at load point G. 

Windage loss (WND) . - The program user may elect to have the program calculate 
an approximate value of windage loss for the main air gap only. Disk losses due to the 
rotor ends are not included. The method used is derived from the results presented in 
references 11 and 12. The windage loss, assuming that the alternator rotor is cylindri- 
cal, is calculated first. In most cases this assumption is valid. The rotors of many 
aircraft salient -pole alternators are made nearly cylindrical by inserting metal pieces 
or coil wedges between the poles at the rotor outside diameter. To account for the effect 
of stator slots and the nonsmoothness of the rotor, the windage loss value is multiplied 
by 1. 5. In calculating the windage loss, the user has the option of assuming either con- 
stant or variable pressure in the air gap cavity as a function of load. 

Efficiency (E) . - At each load efficiency is calculated from 

... . alternator power output x 100 

efficiency = - 

(alternator power output) + (losses) 


Reactances and Time Constants 

In the program the following reactances are calculated: 

(1) Armature winding leakage (XL) 

(2) Direct-axis armature reaction (XD) 

(3) Quadrature-axis armature reaction (XQ) 

(4) Direct-axis synchronous (XA) 

(5) Quadrature -axis synchronous (XB) 

(6) Direct-axis transient (XU) 

(7) Direct-axis subtransient (XX) 

(8) Quadrature-axis subtransient (XY) 

(9) Negative sequence (XN) 

(10) Field leakage (XF) 

(11) Direct -axis damper bar leakage (XDD) 

(12) Quadrature -axis damper bar leakage (XDQ). 

Also, the following time constants are determined: 

(1) Open-circuit time constant (field only) (TC) 

(2) Armature time constant (TA) 
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(3) Transient time constant (TTC). 

The equations used to calculate both the reactances and time constants are given in ref- 
erence 4. The majority of these equations were derived from the theory presented in 
reference 13. 


CONCLUDING REMARKS 

This report describes a digital computer program for the analysis of salient, wound 
pole alternators. The program accepts as input a complete electromagnetic design; 
from this, it calculates the open-circuit saturation curve, field-current requirements at 
rated voltage for various loads, losses, and efficiencies, several reactances and time 
constants, and weights of electromagnetic components. 

The program is based on previous programs written for the Lundell and homopolar 
inductor alternators. Many of the calculation procedures used in this program are based 
on techniques developed for these two alternators. However, the method used in mag- 
netic calculations for the salient, wound pole alternator had to be derived especially for 
this machine. This method, including leakage flux calculations, is described in detail 
in the report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 9, 1973, 

503-35. 
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APPENDIX A 


HOWTO USE COMPUTER PROGRAM 

Input Data Requirements 

To use this computer program for the analysis of salient, wound-pole alternators, 
the complete electromagnetic design of the alternator must be known. This includes 
physical dimensions, armature and field winding parameters, and the magnetic charac- 
teristics of the materials to be used in the stator and rotor. The design information 
must then be transferred to data cards for use with the program. A typical set of data 
cards is shown in figure 7. It consists of two material decks followed by any number of 
alternator design decks. The material decks must be in the order shown in the figure, 
that is, stator material and rotor material. There must be exactly two material decks 
in each data deck even if both materials are identical. The last two cards of the alter- 
nator design deck (the air gap gas name and $WIND) are needed only if windage loss is 
to be calculated by the program. 

If more than one alternator design deck is included in the data deck, the program 
will treat each design deck independently. Each will result in a separate alternator 
analysis complete with an individual output record. However, the same two material 
decks will be assumed to apply to each alternator design deck. 


Preparation of Material Decks 

A material deck consists of five cards. The first card contains the material name. 
This serve two functions: it identifies the material deck, and it is read by the computer 
and stored for later printout on the output record. The remaining four cards contain the 
coordinate values of 14 arbitrary data points located on the magnetization curve of the 
material specified on the first card. 

During program execution, the original magnetization curve is approximately re- 
constructed by interpolation between points. The interpolation assumes a straight line 
on semilog paper between data points. 

Table II lists those variables associated with the two material decks. It also states 
the FORMAT used on each card and gives a brief description of the data appearing on 
each data card. 

To illustrate preparation of a material deck, SAE 4340 steel will be used as an 
example. The first card of this material deck will appear as shown in figure 8. The 
material name should start in column 1 and may extend up to column 36. 
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To prepare the remaining four cards of the material deck, the magnetization curve 
of the material is needed. The magnetization curve for SAE 4340 steel is shown in fig- 
ure 9. The units must be kilolines per square inch for the magnetic flux density and 
ampere-turns per inch for the magnetizing force. Fourteen points on the curve must 
then be chosen. In the figure, 13 points are indicated by data symbols; the 14th point is 
off the graph. These points are listed in the table insert. Careful attention must be paid 
to the sequence in which the numbers are punched onto data cards. The first number 

must be the maximum allowable flux density for the material (B ). During program 

max 

execution, if a flux density exceeds B_ . the alternator is assumed to have saturated, 

max 

and magnetic calculations will be terminated. In the example, this value is 122. 5 kilo- 
lines per square inch. The maximum allowable flux density is followed in ascending 
order, by alternate values of magnetic flux density and magnetizing force. Again, in the 
example, with reference to the table insert, the values appear in the following sequence 
on the data cards: 122.5, 0, 0.01, 6.45, 12, 12.9, 24, . . . , 116, 280, 122.5, 454. 

The complete material deck for SAE 4340 steel is shown in figure 8. 


Preparation of Alternator Design Deck 

The alternator design deck contains all the dimensions, the geometric configuration 
(in numerical code), and the winding parameters needed for an electromagnetic analysis 
of the alternator design. Unlike the material decks, which are read according to a FOR- 
MAT statement, the alternator design decks (except for the air-gap gas name card) are 
read with a READ statement referencing a NAMELIST name. For each NAMELIST name 
one or more data cards are required to numerically define the variables included in that 
NAMELIST name. In all there are nine NAMELIST names. Each name is suggestive of 
the type of variables included in its list. Detailed information about each NAMELIST 
name is provided in table III, which lists all variables used with the alternator design 
deck. All variables belonging to the same NAMELIST name are grouped together. The 
NAMELIST names are arranged in the order in which the data cards must appear in the 
data deck. Units are given, where applicable, and each variable is classified as man- 
datory (M), conditional (C), or optional (O). A mandatory classification indicates that 
the variable must be read in. The conditional classification indicates that, for some 
alternator designs, the variable is required and that, for others, it may be omitted from 
the program input. Variables identified as optional are read in at the discretion of the 
program user. In each case where an optional variable is omitted from the program in- 
put, an assumption regarding that variable is made within the program. The remarks 
column of the table supplies specific explanations in regard to all conditional or optional 
variables. 
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The air-gap gas name card uses an input format of 6A6, which is identical to that 
used for the stator and rotor material names (see table II). The gas name (FLUID) 
should start in column 1 and may extend up to column 36. 

To further clarify the definition of some variables, figures 10 to 14 are given. 
These figures are referenced in table III where applicable. 

Preparation of an alternator design deck will be illustrated with the construction of 
a typical data card for the NAMELIST name DAMPER. The data used are the same as 
that used to obtain the sample output presented in appendix B. Figure 15(a) gives all 
the pertinent design data. Figure 15(b) shows how the design data are related to the 
variables of NAMELIST name DAMPER and how these data are transferred to the data 
card $DAMPER. 

Data cards for the remaining NAMELIST names are prepared in a similar manner. 
To illustrate the results, a complete data deck listing for the sample alternator design 
is shown in figure 16. The output that resulted from this data deck is presented in ap- 
pendix B. 
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APPENDIX B 


TYPICAL OUTPUT LISTING 

The following is the output listing of the salient, wound pole computer program re- 
sulting from the input data shown in figure 16. While this output listing is typical, the 
output format may vary somewhat, depending for example, on the type of stator slot 
configuration used or on whether windage loss was computed by the program. 

Note that on page 21 under "No-Load Saturation Data, " the last column contains all 
zeros. In general, when one or more columns contain all zeros, it indicates that the 
flux density calculated for some section of the magnetic circuit was greater than the 
maximum flux density given on the appropriate material data deck; that is, the alterna- 
tor has saturated thereby causing the calculations to be terminated. When this occurs a 
message is printed on the affected page of the output identifying which section of the al- 
ternator saturated. 


444 WOUND-POLE, SALIENT-POLE ALTERNATOR 444 


ALTERNATOR RATING 


ALTERNATOR KVA 200.0 

LINE-LINE VOLTAGE 294. 

LINE-NEUT. VOLTAGE 170. 

PHASE CURRENT 392.16 

POWER FACTOR 0.75 

PHASES 3 

FREQUENCY 500. 

PULES 12 

RPM 5000.0 


STATOR SLOTS 


TYPE-PARTI ALLY CLOSED 

-BU- 

* 4 

HO * 4 

* * 

BO 0.070 INCHES * * 

OS 0.142 HT * 4 

HO 0.030 4 4 

HX 0.345 4 4 

HT 0.055 4 4 

HW 0.015 HW 4 4 

HS 0.455 4 44444444 4 HS 

4 4 4 4 

NO. OF SLOTS 180 4 4 4 4 

4 44444444 4 

SLOT PITCH 0.218 INCHES HX 4 4 

4 44444444 4 

SLOT PITCH 44 44 

AT 1/3 DIST. 0.223 INCHES 44 44 

4 44444444 4 

4 4 

444444444444 44 

1 1 

1 BS 1 

1 1 
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AIR GAP 


MINIMUM AIR GAP 
MAXIMUM AIR GAP 
EFFECTIVE AIR GAP 


0.070 INCHES 

0.070 

0.076 


CARTER COEFFICIENT 

STATOR 1.066 

ROTOR 1.013 


ARMATURE WINDING ( Y-CONNECTED , FORM WOUND) 


STRAND DIMENSIONS 0.1120 

DISTANCE BTWN CL CF STRANDS (RADIAL) 0.1600 

STRANOS/CONDUCTOR IN RADIAL DIR. 1. 

TOTAL STRANDS/CONDUCTOR 1. 

CONDUCTOR AREA 0.0166 

CURRENT DENSITY AT FULL LOAD 7879.38 

COIL EXTENSION BEYOND CORE 0.250 

MEAN LENGTH OF 1/2 TURN 10.177 

END TURN LENGTH 4.977 

END TURN EXTENSION 1.702 

STATOR SLOT SKEW 0. 

RESISTIVITY AT 20 DEG. C 0.6940 

STATOR RESISTANCE AT 25. DEG. C 0.0058 

NO. OF EFFECTIVE SERIES TURNS 18.20 

TOTAL EFFECTIVE CONDUCTORS 114.13 

WINDING -PITCH 0.8000 

SLOTS SPANNED 12. 

SLOTS PER POLE PER PHASE 5.00 

CONDUCTQRS/SLOT 2. 

NO. OF PARALLEL CIRCUITS 3. 

PHASE BELT ANGLE 60. 

SKEW FACTOR 1.000 

DISTRIBUTION FACTOR 0.957 

PITCH FACTOR 0.951 


FIELO WINDING 

CONDUCTOR DIMENSIONS 0.1900 

CONDUCTOR AREA 0.0047 

NO. OF TURNS (PER POLE) 44. 

MEAN LENGTH OF TURN 15.233 

RESISTIVITY AT 20 DEG. C 0.6940 

FIELD RESISTANCE AT 25. DEG. C 1.1982 

(COILS IN SERIES) 

COIL HEIGHT 0.750 

COIL WIDTH 0.420 


X 0.1500 INCHES 

SQ- IN. 

AMP/ SQ- IN. 

INCHES 

MICRO OHM INCHES 
OHMS 


DEGREES 


X 0.0250 INCHES 
SQ-IN. 

INCHES 

MICRO OHM INCHES 
OHMS 
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CONSTANTS 


Cl, FUNDAMENTAL/ NAX„ OF FIELD FLUX 1,096 
CP, POLE CONSTANT 0.697 
CM, DEMAGNETIZATION FACTOR 0.855 
CO, CROSS MAGNETIZATION FACTOR 0.458 
Dl, POLE FACE LOSS FACTOR 1.170 


STATOR 


STATOR INSIDE DIAMETER 
STATOR OUTSIDE DIAMETER 
OVERALL CORE LENGTH 
EFFECTIVE CORE LENGTH 
DEPTH BELOW SLOT 


12.50 INCHES 

14.50 
5.20 
4.99 
0.55 


STACKING FACTOR 


0.96 


NO. OF COOLING DUCTS 
WIDTH OF DUCTS 


0 . 

0. INCHES 


CORE LOSS AT 77.4 K ILGL INES/SQ. I N. 11.8 WATTS/LB. 
LAMINATION THICKNESS 0.006 IN. 


ROTOR 


POLE 

BODY 

WIDTH 

1.300 

INCHES 



AXIAL LENGTH 

5.200 




HEIGHT 

0.750 




STACKING FACTOR 

1.000 


POLE 

HEAD 

WIDTH 

2.146 

INCHES 



AXIAL LENGTH 

5.200 




HEIGHT 

0.480 




STACKING FACTOR 

0.900 




LAMINATION THICKNESS 

0.006 

INCHES 

POLE 

EMBRACE 

0.667 


ROTOR 

DIAMETER 

12.360 


PERIPHERAL 

SPEED 

16192. 

FEET/MIN. 

SPEC. 

TANGENTIAL FORCE 

2.077 

LBS/SQ. IN. 

DIAMETER OF ROTOR CORE 

9.900 

INCHES 


INSIDE DIAMETER (OF HOLLOW SHAFT) 8.200 
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DAMPER BARS (ROUNO) 


DAMPER BAR DIAMETER 
SLOT OPENING WIDTH 
SLOT OPENING HEIGHT 
DAMPER BAR LENGTH 
DAMPER BAR PITCH 

NO. OF DAMPER BARS/PCL E 

RESISTIVITY AT 20 DEG. C 

TEMPERATURE (HOT) 


0.100 INCHES 
0.0 30 
0.070 
5.2 00 
0.213 

9 

0.694 MICRO-OHM INCHES 
150.00 DEG.C 


WEIGHTS 

STATOR COND. 19.511 POUNDS 

FIELD COND. 12.264 

STATOR IRON 44.539 

ROTOR 67.545 

TOTAL 

(ELECTROMAGNETIC) 143.859 


PERMEANCES 

WINDING LEAKAGE 

(PER INCH OF CORE LENGTH) 

STATOR SLOT 2.998 LINES/AMPERE TURN 

STATOR END 3.456 

OAMPER-O I RECT (PER INCH) 5.772 

DAMPER-QUADRATURE 4.575 


PI POLE HEAD INNER 
P2 POLE HEAD END 
P3 POLE BODY INNER 
P4 POLE BODY END 
EFFECTIVE POLE LEAKAGE 


REACTANCES 

AMPERE CONDUCTORS/INCH 1139.604 

REACTANCE FACTOR 1.985 

STATOR WINDING LEAKAGE 
ARM. REACTION (DIRECT) 

ARM. REACTION (QUAD. ) 

SYNCHRONOUS (DIRECT) 

SYNCHRONOUS (QUAD.) 

FIELD LEAKAGE 
TRANSIENT 

DAMPER LEAKAGE (DIRECT) 

DAMPER LEAKAGE (QUAD.) 


12.813 PERCENT 
163.503 
79.912 
176.316 
92.725 
28.997 
41.810 
11.459 
9.083 


16.919 

3.165 

16.167 

2.567 

58.900 
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SUBTRANSIENT (DIRECT) 
SUBTRANSIENT (QUAD.) 
NEGATIVE SEQUENCE 

FIELD SELF INDUCTANCE 

□PEN CIRCUIT TIME CONSTANT 
(FIELD ONLY, AT 25. DEG. 
TRANSIENT TIME CONSTANT 
ARMATURE TIME CONSTANT 

(WINDING AT 25. DEG. C) 


24.272 

21.896 

23.084 

0.130 HENRIES 


C) 0.10855 SECONDS 

0.02574 

0.01270 


SHORT CIRCUIT AMPERE-TURNS 1479.691 

SHORT CIRCUIT RATIO 0.642 


WINDAGE AT RATED LOAD 

FLUID PROPERTIES CF AIR 

1 . 210E-05 LBS/SEC-FT 
0.0549 LBS/CU . FT . 
15.00 LBS/SQ. IN 
150.00 DEG. C 
29.90 


VI SCOSITY 
DENSITY 
PRESSURE 
TEMPERATURE 
MOLECULAR WEIGHT 


REYNOLDS NUMBER 7139.84 


WINDAGE LOSS 


150.49 WATTS 


STATOR MATERIAL - VANADIUM PERMENDUR 

ROTOR MATERIAL — SAE 4340 STEEL (ROOM TEMP) 

MAGNETIZATION CHAR ACTER I S T I CS (NO LOAD, RATED VOLTAGE) 


TG, TOTAL USEFUL FLUX 

7234.08 KILOLINES 

FQ, MAIN AIR-GAP FLUX/POLE 

420.44 

PPRL , LEAKAGE FLUX/POLE 

50.26 

FLUX DENSITIES ( K I LOLI NE S/ SQ IN) 


MAIN AIR GAP 

35.426 

STATOR CORE 

77.268 

STATOR TEETH 

103.345 

ROTOR - POLE 

69.630 

CORE 

53.2.47 

AMPERE-TURNS/POLE 

MAIN AIR GAP 

839.225 

STATOR CORE 

5.819 

STATOR TEETH 

1.999 

ROTOR - POLE 

54.565 

CORE 

49.028 

TOTAL 

950.636 
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ALTERNATOR LOAD CHARACTERISTICS (RATED VOLTAGE, 0.75 POWER FACTOR) 


PERCENT KVA 

0 . 

50. 

100. 

150. 

200. 

LEAKAGE FLUX/POLE (KL) 

50.26 

87.46 

130. 56 

175.69 

22? .08 

FLUX DENSITIES (KL/SQ-IN) 

STATOR - CORE 

77.27 

82.95 

88.91 

95.99 

104 .07 

TEETH 

103.35 

110.94 

118.92 

128.39 

139.19 

ROTOR - POLE 

69.63 

79.71 

90.88 

103.26 

1 16 .62 

CORE 

53.25 

60.95 

69. 50 

78.96 

89 .18 

AM PE RE-TURNS /POLE 

RAIN AIR GAP 

839.22 

^ 900 . 9 3 

965.72 

1042.58 

1130.31 

STATOR - CORE 

5.82 

6.05 

6.30 

6.99 

a. os 

TEETH 

2.00 

2.50 

3.41 

5.50 

14 .37 

ROTOR - POLE 

54.57 

60.01 

76.33 

120.86 

251 .99 

CORE 

49.03 

52.41 

56.49 

62.80 

82.30 

DEMAGNET 1 7 I NG 

0 . 

568.70 

1234.05 

1919.85 

2607.52 

TOTAL AMPERE TURNS 

950.64 

1590.59 

2342.30 

3158.57 

4094 .58 

FIELO AMPERES (PER COIL) 

21.61 

36.15 

53.23 

71.79 

93.06 

CURRENT DENS. (FIELO) 

4548.50 

7610.49 

11207.20 

15112.77 

19591 .28 

FIELD VOLTS (PER COIL) 

2.78 

5.15 

8.90 

14.71 

24 .37 

TEMPERATURES (DEG.C) 

FIELO 

100.00 

135.49 

200.00 

297. 98 

446.13 

ARMATURE 

100.00 

125.00 

200.00 

325.00 

500 .00 

RESISTANCES (OHMS) 

FIELD (PER CUIL) 

0.13 

0.14 

0.17 

0.20 

0 .26 

ARMATURF 

0.0075 

0.0080 

0.0097 

0.0125 

0.0 164 

ALTERNATOR LOSSES (WATTS) 

FIELD 

721.01 

2232.68 

5686.00 

12671.38 

27219.13 

WINDAGE 

0 . 

101.86 

150.49 

195. 49 

241 .81 

STATOR TOOTH 

595.06 

685.78 

787.97 

918.38 

1079.46 

STATOR CURE 

1192.10 

1373.84 

1578.55 

1839.82 

2162.50 

POLE FACE 

114.97 

133.10 

187.52 

278.20 

405.16 

DAMPER 

0.02 

0.02 

0.03 

0.05 

0 .07 

STATOR cupper 

0 . 

924.89 

4471.48 

12955.52 

30236.58 

EDDY 

0 . 

168.58 

557.90 

974.81 

1320.08 

RISC. LOAD 

0 . 

1000.00 

2000.00 

3000.00 

4000 .00 

TOTAL 

2623.14 

6620.74 

15419.92 

32833.60 

66664 .71 

ALTERNATOR OUTPUT (KVA) 

0 . 

100.00 

200.00 

300.00 

400 .00 

ALTERNATCR OUTPUT (KW) 

0 . 

75.00 

150.00 

225.00 

300 .00 

EFFICIENCY (PER CENT) 

ELECTRO-MAGNETIC 

0 . 

92.00 

90.76 

87.33 

81 .87 

OVER-ALL 

0 . 

91.89 

90.68 

87.27 

81.82 



ALTERNATOR LOAO CHARACTERISTICS (RATED VOLTAGE, I. 00 POWER FACTOR) 


PERCENT KVA 

0. 

50. 

100. 

150. 

200 . 

LEAKAGE FLUX/POLE (KL) 

50.26 

70.37 

105.19 

147.71 

192 .90 

FLUX DENSITIES (KL/SQ-IN) 

STATOR - CORE 

77.27 

82.25 

78.22 

77.94 

81 .54 

TEETH 

103.35 

110.01 

104.61 

104.25 

109.06 

KOTOR - POLE 

69.63 

76.62 

78.52 

84.59 

94.17 

CORE 

53.25 

58.59 

60.04 

64.68 

72 .01 

AMPERE-TURNS/POLE 

PAIN AIR GAP 

839.22 

893.37 

849.53 

846.54 

885.59 

STATOR - CORE 

5.82 

6.02 

5.86 

5.85 

5 .99 

TEETH 

2.00 

2.42 

2.07 

2.05 

2 .35 

ROTOR - POLE 

54.57 

57.87 

58.94 

64. 39 

84.36 

CORE 

49.03 

51.53 

52.07 

53.85 

57.91 

DEMAGNETIZING 

0. 

286.40 

921.84 

1646.58 

2373.94 

TOTAL AMPERE TURNS 

950.64 

1297.60 

1890.30 

2619.26 

3410 .14 

FIELD AMPERES (PER COIL) 

21.61 

29.49 

42.96 

59.53 

77 .50 

CURRENT DENS. (FIELD) 

4548.50 

6208.62 

9044.51 

12532.32 

16316.45 

FIELD VOLTS (PER COIL) 

2.78 

3.99 

6.49 

10.64 

16.96 

TEMPERATURES ( DEG.C ) 

FIELD 

100.00 

117.02 

158.25 

229.98 

334.04 

ARMATURE 

100.00 

125.00 

200.00 

325.00 

500.00 

RESISTANCES (OHMS) 

FIELD (PER COIL) 

0.13 

0.14 

O'. 15 

0. 18 

0 .22 

ARMATURE 

0.0075 

0.0080 

0.0097 

0.0125 

0.0164 

ALTERNATOR LOSSES (WATTS) 

FIELD 

721.01 

1411. 74 

3347.40 

7600.83 

15770 .54 

WINDAGE 

0. 

101.86 

150.49 

195.49 

241 .81 

STATOR TOOTH 

595.06 

674. 32 

609.76 

605.49 

662 .64 

STATOR CURE 

1192.10 

1350.89 

1221.55 

1212.99 

1327.48 

POLE FACE 

114.97 

133.10 

187.52 

278. 20 

405.16 

DAMPER 

0.02 

0.02 

0.03 

0.05 

0 .07 

STATOR CUPPER 

0. 

924.89 

4471.48 

12955.52 

30236 .58 

EDDY 

0. 

168.58 

557.90 

974.81 

1320.08 

PISC. LOAD 

0. 

1000.00 

2000.00 

3000.00 

4000 .00 

TOTAL 

2623.14 

5765.38 

12546.11 

26823.32 

53964 .28 

ALTERNATOR OUTPUT (KVA) 

0. 

100.00 

200.00 

300.00 

400.00 

ALTERNATOR OUTPUT (KW) 

0 . 

100.00 

200.00 

300.00 

400 .00 

EFFICIENCY (PER CENT) 

ELECTRO-MAGNETIC 

0. 

94.64 

94. 16 

91.85 

88.16 

CVER-ALL 

0. 

94.55 

94. 10 

91.79 

88.11 



NO- LOAD SATURATION DATA 


VOLTAGE 


PERCENT 

70. OC 

85.00 

100.00 

115.00 

130.00 

145.00 

147.00 

149.00 

L I NE-NEUTRAL 

119.00 

144.50 

170.00 

195.50 

221.00 

246.50 

249.90 

253.30 

L I NE-L INE 

206. 1 1 

250.28 

294.45 

338.6? 

382.78 

426.95 

432.84 

438.73 

FIELD CURRENT 

15.50 

18.54 

21.61 

24.72 

28.17 

38.09 

44.23 

45.46 

FLUX PER POLE 

294. 3 1 

357.37 

420.44 

483.50 

546. 57 

609.64 

618.04 

626.45 

LEAKAGE FLUX 

35.30 

42.77 

50.26 

57.81 

65.63 

88.30 

102.76 

128.81 

FLLX DENSITIES 
STATOR 
CORE 

54. 09 

65.68 

77.27' 

88.86 

100.45 

112.04 

11 3. 58 

115.13 

TEETH 

72. 34 

87.84 

103.35 

118.85 

134. 35 

149. 85 

151.92 

153.98 

ROTOR 

POLE 

43. 76 

59.19 

69.63 

80.08 

90.56 

103.24 

106.63 

111.73 

CORE 

37.29 

45.27 

53.25 

61.24 

69.25 

78.95 

81.54 

85.44 

AMPERE-TURNS/FOLE 
MAIN AIR-GAP 

587. 46 

713.34 

839.22 

965.11 

1090.99 

1216.88 

1233.66 

1250.44 

STATOR 

CORE 

5.13 

5.47 

5.82 

6.30 

7.52 

10.34 

10.96 

11.62 

TEETH 

1.43 

1.58 

2.00 

3.40 

8.44 

263.48 

490.96 

914.87 

ROTOR 

POLE 

45.71 

50.07 

54.57 

60.51 

76.12 

122.37 

145.02 

190.95 

CORE 

42.07 

45.42 

49.03 

52.51 

56. 35 

62. 79 

65. 73 

72.37 

TOTAL 

681. 80 

815.87 

950.64 

1087.83 

1239.43 

1675.85 

1946.33 

2440.26 


0. 0 

0 . 0 

0 . 0 

0 . 0 


to 


*** STATOR TOOTH IS SATURATED *** 
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APPENDIX C 


COMPLETE FORTRAN LISTINGS AND FLOW CHARTS OF SALIENT, 

WOUND POLE ALTERNATOR ANALYSIS PROGRAM 


The complete FORTRAN listings of the main program and the four subroutines, 
which together constitute the salient, wound pole alternator analysis program, are con- 
tained herein. The main program is SALENT and the four subroutines are, in the order 
given, SUBSAL, OUTPUT, MAGNET, and WINDGE. Each program listing, is followed 
by its flow chart. 


C 

C 

C 


C 

c 

c 


c 

c 


c 


MAIN PROGRAM SALENT 

FOR USE WITH WOUND POLE, SALIENT POLE ALTERNATORS 

COMMON A, AA, AB, AC, AC ORE, AI , ALPHAE, ALPHAR, ALPHAS, APOLE, APOLHO, ARCOR 
1E»AS,AT00TH,B,B1,B2, B3 , BCL , BCO I L , BG, BK , BN, BO , BPL , BRC , BS, BTL , BV, C , C 
21,CC,CCR,CE,CF,CL,CM,CQ,CP,Dl,DCORE,DD,DF,DI»DISH,DR, DRCORE , DSH, DU 
3,DW,DWl»EC,EE,EL,ENDEXr,EP,EW,F»FCL,FE»FFL,FGLl,FGML,FGXL,FH,FKl,F 
APL,FQ,FRC,FS,FTL,G,GA,GC,GE,GP,H,HC,HCOIL,HD,HM,HO,HS,HT,HV,HW,HX, 
5HV , IBN,IPN»IPX,IQQ,I WF,I ZZ , K SAT , LTR, L TR 1 , L TS , P 1 , P 2 , P 3 , P A, P BA , PBH , P 
6BL , PBW , PC , PDD , PDQ, PE , PEFF, PF , PFC , PHH , PHL , PHW, P I , PN , PPL , PPRL , PT , P TC 
7H,PX,QN,GQ»RC,RD»RE,RF,RG1,RK,RK1,RPM,RR»RS,RT,RY,S,SB,SC»$D,SF,SH 
8, SI , SIGMA, SK,SN,SN1, S S , STAT ET , T 1 , T 1 1 , T 2 , T 22 , T3 , T33 , T A, TB , TC , TF , TG , 
97 S , TST ,TT ,TTC»VA,VMI N,VR,WC,WI , WL , WO , WROTOR, WT OT AL , X A , XB , XD , XDD, XD 
$Q, XF , XL » XN » XQ, XR , XU, X X ,X Y, YY , ZG, ZZ 

INTEGER ZZ 

REAL LTS , LTR, LTR 1 

DIMENSION FQLL(IO), QPERV(IO), QVLL(IO), QVLN(IO), FGLLl(lO), 8PLL 
1(10), FELL ( 10) , BCLL(IO), BTLL(IO), FPLL(IO), FTLL(IO), FCLL(IO), 
2FI (10) , BRCORL ( 10 ) , FRCORL(IO), PPRLL(IO), FDMLLI5), TTB ( 5 ) , T TA ( 5 
3), RRA ( 5 ) , RRB ( 5 ) » EZ(5), WMIS{5), CDD(5>, EF(5), PR(5), ST(5), WQ 
AL ( 5 ) , PP ( 5 > , PS ( 5 ) , E X ( 5 ) , SP(5), AKVA(5), WA(5), El(5), E(5), DL ( 
55), WND ( 5 ) , A I ( 60 ) , G(5), RM AT ( 6 ) , SMAT(6), PFC(5), CMPNT(16) 


DATA CMPNT ( 1 ) /2AH *** ROTOR POLE/ ,CMPNT ( 5 ) /2AH * 
1#* ROTOR CORE/ , CMP NT (9)/2AH *** STATOR TOOTH/ , CMPNT ( 1 3 ) /2 AH 
2 *** STATOR CCRE/ 


READ (5,1) 
READ (5,2) 
READ (5,1) 
READ (5,2) 


SMAT 

( AI ( I ) , 1=1 ,29 ) 
RMAT 

( AI ( I ), 1=31,59) 
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1 FORMAT (6A6) 

FORMAT (8F10.1) 

CALL SUBSAL 
CALL OUTPUT 

COMPUTE TOOTH WIDTH AT 1/3 DISTANCE FROM NARROWEST SECTION 

IF (ZZ-3) 4,5,6 
SM=TT-BS 
GO TO 8 

SM= ( 3. 141 6* ( D I +2. *HS ) /QQ)-B3 
GO TO 8 

IF (ZZ-4) 5,7,4 
SM = T T- o 94*BS 
CONTINUE 

AREAS AND LENGTHS FOR MAGNETIC CALCULATIONS 

ATOOTH=QQ*SS*SM*PE/PX 
ACORE=SS* ( I DU- ( D I +2. *HS ) )/2. ) 

CC0RE=( <2.*( DH-2.*HS ) +DU )/ 3. >*3.1416/(2.0* PX) 
AP0LE=PBL*P8W*RK 

APOLHD=0.25*RK1*(PHL+PBL)*(PHW+PBW ) 

ARCGRE=PBL*( ( DSH-DI SH ) /2 .0 )*RK 
DRCORE=( (2. 0*DSH + DISH)/3.0)*3. 14159/(2. 0*PX) 

INITIALIZE SUBSCRIPTED VARIABLES USED IN LOAD CHARACT. CALCS 

DO 9 J=1 , 230 
FGLL1( J)=0. 


NO-LOAD, RATED VOLTAGE MAGNETIZATION CHARACTERISTICS 


FGXL=0. 

PPL-FQ 
FGL 1=FH 
CALL MAGNET 
J=1 

FPLL ( J I =F PL 
FTLL ( J ) =F TL 
FCLL ( J )=FCL 
FRCORL ( J ) =FRC 
FFLL ( J )=FFL 
FGLL 1 ( J ) =FGL1 
FDMLL ( J) =FGXL 
BPLLt J)=BPL 
BTLL ( J ) = BTL 
BCLL ( J ) =BCL 
BRCORL ( J ) = BRC 
PPRLL( J)=PPRL 

SHORT CIRCUIT RATIO AND SHORT CIRCUIT AMPERE-TURNS CALCS 

FSC=XA*FH*0.01 

SCR = FFL/F SC 

WRITE (6,10) F SC , SCR 


SALENT 
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10 FORMAT (IHK,9X,27H SHCRT CIRCUIT AMPERE-TURNS , F 16. 3/ 10X , 20H SHORT 
1CIRCUIT RATIO, F2 3. 3) 

C 

C WINDAGE LOSS AT RATED LOAD 

C 

IF ( IWF.NE.l) GO TO 11 
KW I ND= 1 

CALL WINDGE ( G , RPM, DR , CL , KW I NO , GC , M, WF TTL ) 

C 

11 WRITE (6,12) SMAT 

12 FORMAT ( 1H1 , 18H STATOR MATERIAL -,1H ,6A6) 

WRITE (6,13) RMAT 

13 FORMAT (1HK.18H ROTOR MATERIAL — ,1H ,6A6) 

WRITE (6,14) TG.FQ.PPRL 

14 FORMAT ( 1HL, 30H MAGNETIZATION CHARACTER I ST I CS , 25H (NO LOAD, RATED 
1VOLTAGE)//10X,22H TG , TOTAL USEFUL FLUX , F20 . 2 , 10H K I LOL I NE S/ 10X, 27 
2H FQ, MAIN AIR-GAP F LU X/ POLE , F 1 5 . 2/ 10X , 24H PPRL, LEAKAGE FLUX/POLE 
3,F18.2//10X,33H FLUX DENSITIES ( K I LOL I NES/SO IN)) 

WRITE (6,15) BG.BCL, 8TL, BPL, BRC 

15 FORMAT (13X.13H MAIN AIR G AP , F27 . 3/ 1 3X , 1 2H STATOR CORE , F28 . 3/ 1 3X , 1 
13H STATOR TEETH, F27. 3/13X, 13H ROTOR - POLE , F 27 . 3/ 2 1 X 5H C0RE,F27.3) 

WRITE (6,16) 

16 FORMAT ( 1HK, 9X , 1 8H AMPERE-TURNS/POLE) 

WRITE (6,15) FGL1,FCL,FTL,FPL,FRC 
WRITE (6,17) FFL 

17 FORMAT (1H ,12X,6H TOT AL ,F 34 . 3/// ) 

IF (KSAT.EQ. 10) GO TO 19 

J=4*KS AT 
I = J — 3 

WRITE (6,18) ( CMPNT ( K ) , K = I , J ) 

18 FORMAT ( 1HK,20X,4A6,40H SATURATED AT NO-LOAD, RATED VOLTAGE ***) 

GO TO 3 

C 

C NO-LOAD POLE-FACE LOSS CALCULATION 

C 

19 GT=RO/GC 
AA=l.75/(GT**1.35)+0.8 
GF=AA*PI*SC/(C*FH) 

D2=BG**2. 5*0.000061 

D3= ( 0. 01 67*QQ*RPM)** 1.65*1 . 5147E-5 
IF (TS-0.9) 20,20,21 

20 D4=TS**1. 285*0. 81 
GO TO 24 

21 IF (TS-1.5) 22,22,23 

22 D4=TS**1. 145*0.79 
GO TO 24 

23 D4=TS**0.79*0.92 

24 IF (GT-1.7) 25,25,26 

25 D5=GT**2.31*0.3 
GO TO 31 

26 IF (GT-3.0) 27,27,28 

27 D5=GT**2. 0*0.35 
GO TO 31 

28 IF (GT-5.0) 29,29,30 

29 D5=GT**1.4*0.676 
GO TO 31 

30 D5=GT**0. 965*1. 38 

31 D6=10.0**(0.932*C1-1 .606 ) 

WN=D1*D2*D3*D4*D5*D6*GA 
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c 

c HOT AND COLD DAMPER BAR LOSS CALCULATIONS 

C 

IF (BN) 32,32*33 

32 WD=0«0 
WU=G.O 
GO TO 55 

33 AA=WO/GE 
VT = 0 

IF (AA) 34,37,34 

34 IF ( AA— 0. 65 ) 35,37,36 

35 VT=ALOG.(10.*AA)* (-0.242 1+0.59 
GO TO 37 

36 VT=0.327-(AA *0.266) 

37 CONTINUE 
FS1=2.0*QN*PN*F 
FS2=2.0*FS1 
M=0 

RM=RE*<1.0+ALPHAE*(T 33-20. ) ) 

GO TO 39 

38 RM=RE*(1.0+ALPHAE*(T3-20.) ) 

39 B ARD=DD 

IF (DD.EQ.O.) BARD=SQRT( 4. 0*H*B/3. 1416) 
AA=(FS1/RM)**0.5*BARD*0. 32 
AB= ( FS2/RM)**0.5*BARD*0. 32 
IF (AA-2.5) 40,40,41 

40 Vi=1.0-0. 15*AA+0.3*AA*AA 
GO TO 42 

41 V1 = A A 

42 IF (AB-2.5) 43,43,44 

43 V2 = l .0-0. 15*AB+0.3*AB*A8 
GO TO 45 

44 V2=AB 

45 IF (H.EQ.O. ) GO TO 46 
IF (H.EQ.B) GO TO 46 
VC=H/(3.0*B*V1) 

GO TO 47 

46 VC=0.75/V1 

47 VS=HD/WO+VT+VC 
VG=TB/GE 

Ql=I.0-( 1.0/ ( ( (BO*0.5/GC)**2.0+1.0)**0.5) ) 
GZ=BO/TS 

Q2=1.05*SIN(QZ*2.844) 

IF (QZ-0.37) 48,48,49 

48 Q3=0.46 
GO TO 50 

49 Q3=0.23*SIN( 10.46*QZ-2. 1 )+0. 23 

50 Q4 = SIN(6.283*TB/TS-1 .571 ) + l,0 
Q5=SIN(12.566*T8/TS-1.571)+1 .0 
IF (H) 52,51,52 

51 AB=0.785*DD*DD 
GO TO 53 

52 AB=H*B 

53 W2=PX*BN*SB*RM*1.246/( AB*1000. ) 
W3=(Q2/(2.0*VS+(VG/Q4) ) )**2.0*V1 
W5=(Q3/(2.0*VS+(VG/Q5) ) )**2.0*V2 
WD= ( TS*BG*Q1*CC i **2. 0*W2*( W3+W5 ) 

M = M+1 

IF ( M- 1 ) 55,54,55 
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WU='*D 
GO TO 38 

55 CONTINUE 
C 

C CALCULATE NO-LOAD » RAT ED VOLTAGE TOOTH AND CORE LOSS 

C 

WT=( SM)*QQ*SS*HS*0.849*( BTLLf Ll/BK) **2.0*WL 
WQ= ( DU-HC )*2.67*HC*SS*(BCLL( I) /BK ) 4*2 . 0*WL 
C 

C ARRANGING LOAD POINTS IN INCREASING ORDER 

C 

CO 57 J=l,4 
IA=5-J 

DO 57 1=1, IA 

IF ( G ( D.GT.GI 1 + 1) ) GO TO 56 
GO TO 57 

56 FOL=G( I ) 

G ( I ) =G ( I + 1 ) 

G ( I + 1 ) =POL 

57 CONTINUE 
G ( 1 ) =0 . 

C 

MM = 5 

DO 58 1=2,5 

58 IF <G(I). GE. 1.0. AND. G(I-1).LT. 0.999) MM=I 
C 

C DEFINE NPF 

C 

DO 59 1=1,5 

59 IF (PFC(I).LT. 0.001) GO TO 60 

1 = 1 + 1 

60 NPF= 1-1 
C 

C — 

C CALCULATE ALTERNATOR LOAD CHARACTERISTICS 

C 

1=0 

61 1=1+1 

PF = PFC I I ) 

CK= 1.0 

IF (PF.GE.0.95) CK=1 .10 
AN=ARCOS ( PF ) 

EZ ( 1 ) = 1 .0 
K=0 
I A = 10 
J = 0 
JA = 5 
C 

C ARMATURE TEMPERATURE AND RESISTANCE CALCULATION 

C 

62 J=J+1 

7 T A ( J ) = ( T 1-T1 1 ) *G ( J ) *G ( J ) + T1 1 

RB= ( 1 . OE— 6 ) 4RS* ( 1.0+ALPHAS*( TT A( J ) -20. ) ) 

PRA< J)=RB*RY 
IF ( J.EQ. 1 ) GO TO 62 
C 

C EDDY FACTOR C ALCUL AT I CNS 
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c 

IF (SH) 63,63,64 

63 EZ(J>=1. 

GO TO 65 

64 AA=0.584+( SN*SN- 1 . 0 ) *0.0625*{SD*CL/(SH*HM/2. > > **2 
AB=(SH*SC*F*AC/(BS*R8*1.CE6) )**2.0 

ET=A A* AB* 0.00335*1 .0 
EB=ET-0.00168*AB 
EZI J)=(ET+EB)*0.5 
C 

65 XI = A TANK { ( XB*G ( J ) ) / 100 . ) + S I N( AN ) ) / ( PF+RRA ( J ) * ( P I / EP ) *G ( J ) ) ) 
EB=XI-AN 

EDD=COS(BB)+(RRA(J)*(PI/EP)*G(J)*COS{XI) )+XA*SIN(XI )*G( J )*0.01 
GXX= (EDD-0.93*XD*SIN ( XI ) *G ( J ) *0.01 ) *CK 
FGXL=FGML*G(J)*SIN(XI ) 

PPL=FQ*GXX 

FGL1=FH*GXX 

CALL MAGNET 

FPLL ( J )=FPL 

FTLL ( J )=FTL 

FCLL ( J ) =FCL 

FRCORL ( J ) = FRC 

FFLL ( J )=FFL 

FDMLL ( J ) =FGXL 

BPLL ( J)=BPL 

BTLL ( J )=BTL 

BCLL ( J )=BCL 

BRCORL ( J ) =BRC 

FGLL 1 ( J ) =FGL 1 

PPRL L ( J ) =PPRL 

IF (KSAT.NE. 10) GO TO 66 

IF (K.NE.O) GO TO 67 

IF (J.LT.5) GO TO 62 

GO TO 71 

66 I A=K SAT 

IF ( (G(J)-G(J-l) ).LE. 0.140) GO TO 68 

G ( J ) =G ( J ) -0 . 10 

J=J-1 

K=K+ 1 

GO TO 62 

67 IF (J.EQ.5) GO TO 71 

J = J+1 

68 JA= J— 1 

69 V=220+ J 

DO 70 K=J,M,5 

70 FGLL1 { K) =0. 

J = J+1 

IF (J.LT.6) GO TO 69 
IF (MM.GT.JA) MM= J A 

71 IF (I.EQ.l) FIMM=FFLL(MM)/PT 

VV=3.*PI*EP*PF 

C 

C LOSSES AND EFFICIENCY UNDER LOAD 

C 

M=0 

KWI ND=0 

vw=wu 

72 P = M + 1 
UA=G ( M) 
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FI (M)=FFLL(M)/PT 

CDD(M)=FI(M)/AS 

IF ( MM. NE . 1 ) GO TO 73 

TIB ( M) =T22 

GO TO 74 

73 TTB ( M) = { ( ( T2-T22 )*FI ( M ) )+(T22*FIMM**2-T2*FI ( 1 )**2 ) ) / ( F I MM**2-F I 

1(1) **2 ) 

7 A RRB(M)=(1„0E-6)*RR*< 1 . 0+ AL PH AR* ( TTB ( M ) - 20 . ) )*ZG 
PR ( M ) =F I ( M ) *F I ( M ) *RR B ( M ) 

RRB(M)=RRB(M)/PX 
EF(M)=FI ( M)*RRB(M) 

PS(M)=(3„*(PI*UA)**2 )*RRA(M) 

WQL(M)=WQ#(BCLL(M)/BCLL(1) )**2 
ST ( M)=WT*(8TLL( M)/BTLL( 1 ) ) **2 
WA ( M)=VV*UA/1G00. 

IF ( IWF.NE. 1 ) GO TO 75 
IF ( M.EQ. 1 ) GO TO 75 

CALL WINDGE ( G , RPM , DR , CL , KW I ND , GC , M, WFTTL ) 

WND ( M) =WFTTL 

75 AKVA(M)=WA(M)/PF 

WMI S (M)=AKVA( M)*10.0 
GM=(GF*UA>**2+1.0 
PP ( M ) =GM*WN 
DL ( M )=GM*WW 

EX(M.) = (EZ(M)— 1.0)*PS(M)* (CL/ HM ) 

SP(M) = PP(M)+PR(M)+PS(M)+EX(M)+ST(M)-»-WQL(M)+WMIS(M)+WND(M) 
E(M)=(WA(M)/{WA(M)+(SP{M)-WND<M) )/1000. ) 1*100. 
E1(M)=(WA(M)/(WA(M)+SP(M>/1000. ) )*100. 

IF (M.EQ.l) WW=WD 
IF (M.LT.JA) GO TO 72 

WRITE (6,76) PF, (G( I ) » 1=1,5) , (PPRLL( I ) , 1=1,5 ) 

IF (KSAT.EQ.10) G ( JA ) = G ( JA ) +FLOAT ( K ) *0. 10 
IF (KSAT.NE.10) G ( J A + 1 ) = G( JA +1 ) +FLOAT ( K ) *0 . 1 0 

WRITE (6,77) ( BCLL ( I ) , 1 = 1, 5) , CBTLU I ) , 1 = 1,5) , ( BPLL ( I ) , 1=1,5) , (BRCO 
1 R L ( I ), 1=1,5) 

WRITE (6,78) (FGLL1( I >,1=1,5) 

WRITE (6,77) (FCLL( I ) , 1=1,5) , ( FTLL ( I ) , I = 1 , 5 ) , ( F PLL ( I ) , I = 1 , 5 ) , ( FRCO 
1RL( I ) , 1=1,5 ) 

WRITE (6,79) ( FDMLL( I ) , I =1 , 5 ) , ( FFLL ( I ) , I =1 , 5 ) , ( F I ( I ) , I =1 , 5 ) , ( CDD ( I 

1 ) , 1 = 1,5) , (EF( I ) , 1 = 1,5 ) 

76 FORMAT ( 1H1 , 26X47HALTERNAT0R LOAD CHARACTERISTICS (RATED VOLTAGE, F 

15.2, 14H POWER FACTOR ) /27X, 1 1 ( 6H ) , //7X1 1HPERCENT KVA,12X,2PF1 

27.0, AF 19.0// 7X22HLEAK AGE FLUX/ POLE ( KL ) , 1X0PF 1 9 . 2, AF 19 . 2// 7X25HF LU 
3 X DENSITIES (KL/SQ-IN) ) 

77 FORMAT ( 10X13HST ATOR - CORE , 7X5 F 19 . 2 / 19X 5HTEETH6X 5F 1 9. 2 / 10X , 1 2HRQT 

I OR - POLE, 9X,5F19.2/18X, 4HC0RE , 8X , 5F 19 . 2 ) 

78 FORMAT ( 1H / 7X , 1 7H AMPERE-T URNS/ POLE/ 10X , 1 2HM A I N AIR GAP , 8X , 5F 19. 2 ) 

79 FORMAT ( 1H , 9X 1 3HDEM AGNE TI Z I NG7X5F 19 . 2/7X1 8HT0T AL AMPERE TURNS 5X5F 
119.2//7X24HFIELD AMPERES (PER CO I L ) , F 18 . 2, AF 19 . 2/7X2 1HCURRENT DENS 
2. (FIELD)2X5F19.2/7X22HFIEL0 VOLTS (PER COIL ) , 1X5F19.2 ) 

WRITE (6,80) ( TT8( I ) , 1 = 1,5 ) , (TTA( I ) , 1 = 1,5) , ( RRB ( I ) , 1=1,5) , ( RRA ( I ) , 

II = 1, 5), (PR(I), 1 = 1,5), (WNU( I), 1 = 1,5), (ST(I), 1 = 1,5), (WQL(I)» 1 = 1, 5)»( 
2FP( I ) , 1 = 1 ,5) , ( DL( I ), 1 = 1, 5) , ( PS( I), 1 = 1,5 ), (EX( I ) , 1 = 1,5), ( WM I S ( I ), 1 = 
31,5) ,( SP( I) , 1=1,5) , ( AKVA(I ),I=1,5),(WA( I )»I=1,5)»(E( I), 1=1, 5) 

80 FORMAT ( 1HK , 6X20HTEMPERA TURE S ( DEG . C ) / 10X5HF I ELD 1 5X5F1 9. 2/ 1 0X8HARM 
1 ATURE12X5F19.2/7X18HRESI STANCES ( OHMS ) / 1 OX 16HF I ELD (PER C0IL)4X,5F 
219. 2/1 0X8HARMATURE1 2 X5F 19. A/ /7X25H ALTERNATOR LOSSES ( WATTS )/ 10X5HF 
3 I ELD 15 X5F 19. 2/10X7HW I ND AGE 13X5 FI 9. 2/ 10X1 2HST ATOR TOOTHS X5F 1 9 . 2/ 1 OX 


28 


SALENT 



4 1 1HSTA TOR CGRE9X5F19 . 2/ 1 OX9HPOL E F AC E 1 1 X5F 1 9 . 2/ 1 OX 6HDAMP ER , 1 4X5F 1 9 
5.2/10X13HSTAT0R COPP ER7X5F 19 * 2/ 10X4HF0DY 1 6X 5F 1 9 . 2/ 10X 10HM I SC . LOAD 
610X5F19.2/10X5HTOTAL 1 5X5 FI 9 „ 2/ / 7X2 3HALT E RN ATOR OUTPUT ( K VA ) 5F 1 9. 2/ 
77X22HALTERNATOR OUTPUT ( KW ) 1 X5 F 19 . 2 / /7X , 2 IHEFF I C I ENC Y (PER CENT) /I 
80X , 16HELECTR0- MAGNET I C , 4 X, 5F19.2! 

IF (IWF.EQ.l) WRITE (6,81) (El(l), 1 = 1,5) 

81 FORMAT (10X,8H0VER-ALL»12X,5F19.2) 

KSAT = 1 A 

IF (KSAT.EQ. 10) GO TO 83 
J = 4*KS AT 
M= J- 3 

WRITE (6,82) ( CMPNT ( K ) ,K=M, J ) 

82 FORMAT ( 1HK , 20X , 4A6, 1 3H IS SATURATED, 4H ***) 

83 IF ( I.GE.NPF) GO TO 84 
GO TO 61 

C 

C INITIALIZE VARIABLES USED IN NO-LOAD MAGNETIC CALCS 

C 

84 CONTINUE 

DO 85 J=1 * 160 

85 FQLL ( J ) =0 . 

C 

C 

C CALCULATE NO-LOAD SATURATION DATA 

c 

C 

FGXL=0„ 

I DEL R= 15 
EDD=VMIN 
J=1 

86 PPL=FQ*EDD 
FGL1=FH*EDD 
CALL MAGNET 

IF (KSAT. NE. 10) GO TO 87 

FQLL(J)=PPL 

CPERV( J)=100.*EDD 

GVLL ( J )=EE*EDD 

CVLN( J)=QVLL( J)/SQRT( 3. ) 

FI ( J)=FFL/PT 

FPL L ( J ) =F PL 

FTLL ( J )=FTL 

FCLL ( J ) =FCL 

FRCORL ( J ) =FRC 

FFLL( J)=FFL 

8PL L ( J ) = B PL 

BTLL(J)=BTL 

BCLL ( J ) =BCL 

BRCORL ( J ) =BRC 

FGLL 1 ( J ) =FGL 1 

P PRLL ( J ) = PPRL 

IF (J.EQ.10) GO TO 90 

J = J + 1 

EDD=EDD+FLOAT ( I DELR ) / 100 . 

GO TO 86 

87 EDD = EDD-FLOAT ( I DELR) / 100 . 

IF ( IDELR.GT .5 ) GO TO 88 
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IF ( IDELR 6 GT.a2I GO TO 89 

GO TO 90 

88 IDELR=( IDELR/6 ! *5 
EOD^EDD+FLOAT ( I DELR) /100 . 

GO TO 86 

89 IDELR=( IDELR/3)#2 
EDD=EDD+FLOAT ( I DELR) /100« 

GO TO 86 

90 WRITE (6,91) ( QPERV ( K ) , K = 1 , 1 0 ) , ( QVLN ( K ) , K= 1 , 10 ) » ( Q VLL ( K ) , K = 1 , 1 0 ) , ( 
IF I (K),K=1 ,10) , (FQLL(K),K=1,10> , (PPRLL(K) ,K=1,10) 

91 FORMAT ( 1H1 , 50X23H NO-LOAD SATURATION DAT A/ 51X23H 

1 //2X7HV0LTAGE/5X7HPERCENT6X10F11.2//5X12HLINE-NEUTRAL 1X10F11. 

22/5X9HLINE-LINE4X10F11 . 2 // 2.X 1 3HF I ELD CURRENT3X 10F 1 1 * 2/ /2X , 1 3HFLUX 
3PER POLE* 3X , IOF11.2/2X, 12H LEAKAGE FLUX , 4X1 OF 1 1 . 2 // 2X , 14HFL UX DENSI 
4TIES ) 

WRITE (6,92) (BCLL(K) ,K=1, 10), ( BTLL ( K ) , K=1 , 10 ) , ( BP LL ( K ) , K= 1 , 10 ) , (B 
lRCORL(K) » K = 1 ,10 ) 

92 FORMAT ( 1H , 4X6HST ATOR/7 X, 4HC0RE , 7 X t 10F1 1 . 2/ 7X , 5HT EETH, 6X , 10F 1 1 . 2/ 
15X»5HROTOR/7X» 4HP0LE »7X»10F11.2/7X» 4HC0RE* 7X, 10F11. 2) 

WRITE 16,93) 

93 FORMAT ( 1HK, IX, 17H AMPERE -TURNS/P OLE ) 

WRITE (6,94) (FGLLH K ) »K=1 ,10 ) 

94 FORMAT <1H ,4X,12HMAIN AIR-GAP, F12«2,9F11»2) 

WRITE (6,92) (FCLL(K) ,K=1, 10) , ( FTLL ( K ) , K=1 , 10 ) , ( FPLL ( K ) , K= 1 , 10 ) , ( F 
1RC0RL(K),K=1,10) 

WRITE (6,95) ( FFLL ( K ) , K- 1 , 10 ) 

95 FORMAT ( 1 HK , 4X , 5HT0T AL , 8 X, 10F1 1 . 2 ) 

IF (KSAT.EQ.10) GO TO 3 

J=4*KS AT 
I- J-3 

WRITE (6,82) (CMPNT(K) ,K=I , J J 

GO TO 3 

END 
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SUBROUTINE SUBSAL 

COMMON A,AA,AB, AC, AC ORE, AI ,ALPHAE, AL PHAR , ALPHAS , APOLE , APOLHD > ARCOR 
1E,AS,AT00TH,B,B1,B2, B 3 , BCL » BCO I L, BG, BK, BN, BO » BPL , BRC , BS, BTL , BV , C » C 
21,CC,CCR,CE,CF,CL,CM,CQ,CP,Dl»DCaRE,OD,DF,DI,OISH,DR,DRCORE,DSH,DU 
3,DW,DWl,EC,EE,EL,ENDEXT, EP, EW,F, FCL,FE,FFL ,FGL1, FGML , FGXL , FH, FK 1 , F 
4 PL, FQ, FRC >FS,FTL»G»GA, GC ,GE , GP , H, HC , HCO I L» HD , HM , HO » HS , HT , HV , HW » H X » 
5HY, 18N,IPN, IPX, IGQ. I WF, 1 ZZ , K S AT , L T R, LTR 1 , LTS , P 1 , P 2 , P 3 , P4 , PB A , PBH , P 
6BL,PBW,PC,PDD,PDQ,PE , PEFF, PF , PFC , PHH, PHL , PHW , P I , PN , PPL , PPRL , PT , P TC 
7H,PX,QN,QQ,RC,RD,RE,RF,RG1 , RK , RK1 » RPM, RR , RS , RT , RY , S , SB, SC, SD , SF , SH 
8, SI , SIGMA , SK,SN,SN1, SS.STATET, T1,T 11, T2 , T22, T3 , T33 , TA, TB , TC , TF , T G, 
9TS, TST,TT, TTC,VA,VMI N, VR,WC, Wl , WL , WO, WROTOR , WTOT AL , XA, XB, XD,XDD, XD 
$Q,XF,XL,XN,XG,XR,XU, XX, XY, YY,ZG,ZZ 
C 

INTEGER ZZ 
C 

REAL LT,LTS,LTR,LTR1 
C 

DIMENSION D A ( 8 ) , DX<6>, DY ( 8 ) , DZ(8), G(5), PFC ( 5 ) , AII60) 

C 

NAMELIST /RATING/ VA , EE , EP , F , RPM , I PX, PFC , G , VM I N/ ST ATOR /DI , DU, CL , HV 

1 , BV, SF,LTS,WL,BK/SLOTS/ZZ,BO,B3, BS , HO, HX , HY, HS , HT , I QQ/W I NDNG/RF, SC 

2, BX,SEP,PTCH,EL,YY,C,DW,SN,SN1,DW1,CE,SD,PBA,SK,T1,RS, ALPHAS,T11,T 
3ST/AIRGAP/GC.GP, I WF/ ROTOR/ C 1 , CP , CM , CO, Dl , RK , LTR , RK 1 , LTR 1 , PHL , PE , PH 
4W,WR0T0R,PHH»DISH,DSH,PBL, PBH, PBW/F I ELD/HCOI L , PT , RD, RT , T2 , BCOI L , TF 
5, T22, RR,AL PH AR,F£/ DAMPER/BN, WO, HD, DD, H, B , SB, TB , T 33 , RE , ALPHAE , T3 

C 

DATA DA,DX,DY, DZ/O. 05, 0.072, 0.125,0. 165,0.225,0.438,0.688,1.5,0.00 

10124.0. 00021.0.00021.0.00084.2*0.00189,2*0.000124, 2*0.00084,0.0018 

29.0. 00335.0.00754.0. 0302.3*0.000124,2*0.00335,0.00754,0.0134,0.030 
32/ 

C 

WRITE (6,1) 

1 FORMAT ( 1H1, 14X,43H*** WOUND-POLE, SALIENT-POLE ALTERNATOR ***) 

DO 2 1=1,5 

2 PFC ( I ) =0. 

ALPHAE=0. 00393 
ALPHAR=0. 00393 
ALPHAS=0. 00393 
B = 0. 

BCOI L=0 . 

BV = 0. 

BX = 0 . 

C1=0 
CM = 0 
CP = 0 
CQ = 0 
CW = 0 
D1 = 0 . 

CD = 0 . 

DI SH=0. 

DSH=0. 

DW 1=0 
E E = 0 . 

EL = 0 

ENDE XT =0 . 

EP=0. 

F = 0 . 
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G( 1 ) =0. 

G ( 2 ) =0 . 75 
G<3)=1.00 
G ( A ) = 1 . 25 
G(5)=1.50 
GP =0 . 

H=0. 

HCO l L = 0. 

HV = 0 . 

I PN = 3 
IPX = 0 
IWF = 0 
LTR=0. 

LTR 1 =0 . 

LTS=0. 

PBA=60. 

PBL = 0 
PBW = 0 
PE = 0 . 

PHL=0. 

PHW=0. 

PN = 3. 

PTCH=0. 

RE =0 . 694 
RK = 0. 

RK 1 =0 . 

RPM=0. 

RR=0.694 
RS=0 . 694 
R T = 0 . 

SB = 0 . 

SEP=0. 01 
SF = 0 . 

SK = Q 
SN=1.0 
T3 = 0 . 

T33=20 . 

T B = 0 . 

TF=25. 

TST=25. 

VMIN=0.7 

WOO. 

WR0TOR=0. 

XF = 0. 

YY = 0. 

READ { 5,RAT I NG ) 

READ (5, STATOR) 

READ (5, SLOTS) 

READ (5 t AIRGAP) 

READ (5, ROTOR) 

READ ( 5 , F I ELD ) 

READ (5, DAMPER) 

PF=PFC ( 1 ) 

IF (EP.EQ.O.) EP=EE/1 .732051 
IF (EE.EQ.O.) EE=EP*1. 732051 
IF (DW1.NE.0.) SH=DW 1 

IF ( IPX.EQ.O. AND.RPM.NE.C. ) IPX=(F*120. )/RPM 

PX= IPX 

IF (RPM.EQ.O.. AND.PX.NE.O. ) RPM= { F*1 20. ) /PX 
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IF (F.EQ.O. ) F=PX#RPM/12C. 

HW=H Y-HO— HT 
C Q= I QQ 

IF ( ZZ.NE. 3) GO TO 3 

Bl=(HO+HT-HS!#(6.283185/GQ)+B3 

B2=B1+(6.283185#HW/QQ) 

BS= ( B2+B3 ) /2 . 

IF (BX.EQ.O.) BX=8S- .015 

3 CONTINUE 

IF (BX.EQ.O.) BX=BS- . 0 15 
PI=(VA*1COO. ) / ( EE#SQRT ( 3 . ) ) 

IF ( ZZ.EQ. 1.0R.ZZ.EQ.5) BO=BS 

IZZ=ZZ 

DB=.25 

IF (YY.EQ.O.) YY=FLOAT(IFIX( (QQ/PX)#PTCH+0.45) ) 

IF (PTCH.EQ.O.) PTCH=YY*PX/QQ 
IF (DU.GE.8.) DB=0. 5 
DR = D 1-2. *GC 

IF (DSH.EQ.O.) DSH=DR-2. 0* ( PHH+PBH ) 

IF (HCOIL.EQ.O. ) HCO I L = PBH 
IF ( GP.EQ.O. ) GP=GC 

IF (PE.EQ.O.) PE=(PX/3.1415927)#(ARSIN(PHW/DR)) 

IF (PHW.EQ.O.) PHW=DR*SIN( 3.1415927#PE/PX) 

IF (PBL.EQ.O.) PBL=PHl 
IF ( PHL.EQ.O. ) PHL=P BL 

FE=2.0*PBL+2.0#( PBW-0 .25 )+3. 14159*( 0.45+BCOIL ) 

HC= ( DU-DI — ?.0#HS)#0.5 
Z Y=0.7#HS 
DO 4 1=1,5 

4 IF ( G ( I ) . GT . 9 . ) G(I)=G(I)/100. 

QN=QQ/(PX#PN) 

CS- Y Y/ ( PN#QN ) 

IF (BN.EQ.O. ) GO TO 5 
IF (SB.EQ.O.) SB=PHL 
IF (T3.EQ.0.) T3 = T1 
IF (DD.NE.O.) XX=DD 
IF (H.NE.O.) XX=H 

IF (TB.EQ.O.) TB-3. 14159 27 #( DR-2.0#W0-XX )*P E / ( PX# ( BN+I.O) ) 

CHECK FOR ERROR CONDITIONS 

IF (ABS(PTCH-YY#PX/QQ).GT. 0.0001) GO TO 6 
GO TO 7 

6 WRITE (6,9) PTCH 
PTCH=YY*PX/QQ 

7 IF (CS.GT .l.O.OR.CS.LT .0.5 ) WRITE (6,10) CS 

IF (EP#EE. EQ.O.. OR. ABS(EE/EP-1. 732051). GT. 0.01) WRITE (6,11) 

IF (PX*F*RPM.EQ.0..0R.ABS(F-PX*RPM/120.).GT.0.1) WRITE (6,12) 

IF (HC.LT.ZY) WRITE (6,13) HC.HS 

IF ( ((3.14159#DSH/PX)-PBW).LT.2.0*BC0IL) WRITE (6,14) 

IF (RT.LT. 1.0E-10) GO TO 8 

IF ( ( (HCOIL#BCOIL)/(RT#RO) ).LE.PT) WRITE (6,15) 

GO TO 16 

8 IF (HC OIL# BC OIL/RD#* 2. LE.0.8573*PT) WRITE (6,15) 

9 FORMAT (1H ,F5.3,30H WINDING PITCH IS NOT POSSIBLE) 

10 FORMAT ( 5 X , 2 7H CS (PER UNIT POLE PITCH) = , F7 . 3/ 10X , 31H CS MUST BE 
1BETWEEN 0.5 AND 1.0) 

11 FORMAT (1H , 38H EITHER PHASE OR LINE VOLTAGE IS WRONG) 

12 FORMAT (1H ,44H FREQUENCY, RPM, OR NO. OF POLES IS IN ERROR) 
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13 FORMAT (1H / 5 X54HDEP T H BELOW SLOT IS LESS THAN 70 PERCENT OF SLOT 
1DEPTH/10X,4HDBS=F8.4/10X,4H SD=F8 . 4 I 

14 FORMAT (1H ,57H FIELD COIL IS WIDER THAN AVAILABLE SPACE AT BASE 0 
IF POLE! 

15 FORMAT (1H ,81H FIELD COIL DIMENSIONS ARE TOO SMALL FOR THE SPECIF 
1 1 ED NO. OF TURNS AND WIRE SIZE) 

C 

C DETERMINE ROTOR AND STATOR STACKING FACTORS 

C 

16 M=1 
STFK=SF 
LT=LTS 
GO TO 19 

17 M = 2 

STF K=RK 
LT = LTR 
GO TO 19 

18 M=3 
STFK=RK1 
LT=LTRl 

19 IF (STFK.NE.O.) GOTO (21, 22, 23), M 
IF (LT.EO.O.) GO TO 20 
STFK=1.Q-(12. 5E— 4/LT ) 

GO TO (21, 22, 23), M 

20 STFK=i .0 

GO TO (21, 22, 23), M 

21 SF=STFK 
GO TO 17 

22 RK=STFK 
GO TO 18 

23 RK1 =STFK 
C 

C CALCULATE POLE FACE LOSS FACTOR 

C 

M=0 

IF (Dl.NE.O. ) GO TO 31 
IF (LTR1.NE.0. ) GO TO 24 
M= 1 

IF (RK1.GT. 0.9999) GO TO 30 
LTR 1 = ( 12.5E-4) / ( 1. 0— RK 1 ) 

24 IF (LTR1-0.G45) 25,25,26 

25 01=1.17 
GO TO 31 

26 IF ( LTR1-0 .094 ) 27,27,28 

27 01=1.75 
GO TO 31 

28 IF (LTR1-0.17) 29,29,30 

29 D1 = 3 . 5 
GO TO 31 

30 Dl=7.0 

31 IF (M.EQ.l) LTR1=0. 

IBN=BN+. 1 

C 

C 

SS=SF*(CL-HV*BV) 

SIGMA=(54.E3/DI**2)*( PF/ SS ) * ( VA/RPM ) 

VR=0.262#DR*RPM 

TP=3.14159*DR/PX 
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1S=3.142*DI/QQ 
IF ( ZZ-4 ) 32,33,32 

32 TT=( ,667*HS+DI )*3.142/CQ 
GO TO 34 

33 TT=(DI+2.0*HO+1.333*BS)*3. 1416/QQ 
C 

C CALCULATE CARTER COEFFICIENTS 

C 

34 IF (ZZ.GT.l. AND.ZZ.LT. 5) GO TO 35 
CC=(5.0*GC+BS)*TS/((5.0*GC+BS)*rS-BS*BS) 

GO TO 36 

35 QC=(4.44*GC+0.75*B0)*TS 
CC =QC/ ( QC-RO*BO ) 

36 IF ( IBN.EQ.O) GO TO 37 

QC = ( 4. 44*GC+0. 75*WO) * TB 
CCR=QC/(QC-WO**2) 

GO TO 38 

37 CCR= 1 . 

38 CONTINUE 
C 

C PITCH FACTOR AND SKEW FACTOR CALCULATIONS 

C 

CF=SIN(YY*1.571/(PN*QN) ) 

IF (SK) 39,39,40 

39 F S= 1 . 0 
GO TO 41 

40 FS = ( SK/TP ) *1. 5707 
FS=( l./FS)*(SIN(FS)) 

C 

C CHECK IF WINDING HAS INTEGRAL NO. OF SLOTS PER POLE PER PHASE 

C 

41 D= 1 . 0 

IF ( PBA.GT.61.0) D=2 . 0 
IC=IFIX(C+0.001) 

IZY=IPX*I PN 
IDM = 0 

42 I DM= I DM+ I ZY 

IF (IQQ-IDM) 44,43,42 
C 

C CALCULATE DISTRIBUTION FACTOR FOR INTEGRAL SLOT WINDING 

C 

43 DF= S IN (1. 571 *D/PN) /( QN*D*S IN ( 1.571/1 PN*QN ) ) ) 

IF < UPX/IC)*IC.NE.IPX) WRITE (6,48) IC 

GO TO 50 
C 

C CALCULATE DISTRIBUTION FACTOR FOR FRACTIONAL SLOT WINDING 

C 

44 1 1 QQ = I QQ 
1 = 2 

45 IF ( ( I ZY/ I ) *I.EQ. I ZY. AND.( I IQQ/I >*I.EQ. I IQQ) GO TO 46 
IF (I.GT.IZY) GO TO 47 

1 = 1+1 
GO TO 45 

46 I Z Y= I Z Y/ I 

1 1 G G = I IQQ/I 
GO TO 45 

47 FNQ= 1 1 QQ 

DF=SIN(1. 571*D/PN) /( FNQ*D*S IN( 1.571/ ( FNQ4PN) ) ) 

IF ( ( IZY/3)*3.EQ.IZY) WRITE (6,49) 
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IF ( ( I PX/ I ZY ) *1 ZY.NE . I PX ) WRITE (6,49) 

I ZY= I PX/I ZY 

IF ( (IZY/IC)*IC.NE.IZY) WRITE (6,48) IC 

48 FORMAT (1H ,I2,35H PARALLEL CIRCUITS ARE NOT POSSIBLE! 

49 FORMAT ( 2H ,41H IMPROPER FRACTIONAL-SLOT WINDING IS USED) 

50 EC=QQ*SC*CF*FS/C 
C 

C COMPUTE ARMATURE CONDUCTOR AREA 

C 

IF ( DW1 ) 51,51,52 

51 AC=G.785*DW*DW*SN1 
GO TO 64 

52 Z Y = 0 .0 

DT=AMI N1 ( DW,DW1 ) 

DG= AMAX1 ! DW, DW1 ) 

53 IF (DT-.05) 56,56,54 

54 J A=0 

55 J A = J A+ 1 

IF ( DT-DA ( J A ) ) 57,57,55 

56 0 = 0 

IF (ZY) 63,63,78 

57 IF (DG-0.188) 58,58,59 

58 CY=DX(JA— 1) 

CZ = D X ( JA ) 

GO TO 62 

59 IF (DG-0.75) 60,60,61 

60 CY = D Y ( JA- 1 ) 

CZ=DY( JA) 

GO TO 62 

61 C Y=DZ ( JA-1 ) 

CZ=DZ( JA) 

62 0=CY+(CZ-CY)*(DT-DA( J A-l ) ) / ( DA ( J A ) -DA( J A-l ) ) 

IF (ZY) 63,63,78 

63 AC= ( DT*DG-D)*SN1 
C 

C CALCULATE END EXTENSION LENGTH 

C 

64 IF (EL) 65,65,75 

65 - IF (RF) 66,66,72 

66 IF (PX-2.0) 67,67,68 

67 0=1.3 

GO TO 71 

68 IF (PX-4.0) 69,69,70 

69 U= 1 . 5 

GO TO 71 

70 0=1.7 

71 EL=3.14l6*U*YY*(DI+HS)/QQ+0.5 
GO TO 75 

72 IF (ZZ.EQ.4) GO TO 73 
0=D I +2 . *H Y+HX 

GO TO 74 

73 U=D I +HS 

74 l U= 3 . 1 416*U/ GQ 
SALPHA=(BX+SEP)/UU 
ALPHA=ARSIN(SALPHA) 

CALPHA=COS( ALPHA) 

STRLEN=3. 1 4 1 6*U*YY / ( GC*C ALPHA) 

EL=2.0*CE+STRLEN+1»5708* (0.54HX+DB) 
ENDEXT=STRLEN*SALPHA/2.0+CE+(HX+DB>/2.0 
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HM=CL+EL 


75 
C 

C CALCULATE STATOR RESISTANCE 

C 

A=P(*$C*CF/(C*TS) 

RY=SC*QQ*HM/{ PN*AC*C*C ) 

RGI= ( 1 .E-6) *RS*( 1.0+ ALPHAS *( TST-20. ) >*RY 
S=PI / ( C*AC ) 

C 

C COMPUTE FIELD CONDUCTOR AREA 

C 

IF (RT) 76,76,77 

76 AS=. 7854*RD*RD 
GO TO 79 

77 Z Y=1 .0 

DT=AMI N1 ( RT , RD ) 

DG=AMAX1!RT ,RD) 

GO TO 53 

78 AS=DT*DG-D 
C 

C COMPUTE FIELD RESISTANCE 

C 

79 ZG= ( PT*FE/AS)*PX 

FKI= ( l.E-6)*RR*( 1.0+ ALPHAR*( TF-20. ) ) *ZG 
C 

C NO LOAD MAGNETIC CALCULATIONS 

C 

GA=3.1416*DI*(CL-HV*BV) 

GE=CC*GC*CCR 

AG=6.3B*DI/(PX*GE) 

IF (Cl) 81,80,81 

80 C 1= (. 649* ALOG( PE) +1.359)*! ( GC/GP ) **0 . 352 ) 

81 CW=(0. 707/1. 732)*Cl*DF 
TG=EE/(CW*EC*RPM)*6.0E6 
BG=TG/GA 

FH=BG*GE/0. 00319 
IF (CP) 82,82,83 

82 CP=PE*( ALOGIGC/TP)*. 0378+1. 191)*! ( GC/GP ) **0. 41 ) 

83 FQ= T G*CP / PX 
C 

C DETERMINE DEMAGNETIZING AMPERE TURNS (RATED LOAD) 

C 

IF (CM) 84,84,85 

84 AA=SIN(3. 142*PE) 

AB=SIN(1.571*PE)*4.0 

CM=(3.142*PE+AA)/AB 

85 CONTINUE 

FGML = . 45* EC*PI*CM*OF/PX 
C 

C REACTANCE FACTOR CALCULATION 

C 

XR=.0707*A*DF/(C1*BG ) 

C 

C SPECIFIC ARMATURE SLOT AND END-TURN LEAKAGE PERMEANCES 

C 

FACTOR=YY/(PN*QN) 

IF (PBA.LT .61. ) GO TQ 86 
FF = . 05 * ( 24.*FACT0R-1 . ) 

IF ( FACTOR. GE. 0.667) FF=.75 
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IF (ZZ.EQ.5) F F= 1 . 

GO TO 87 

FF=. 254(6. 4FACTOR-1. ) 

IF (FACT0R.GE.0.667) F F= .2 5* { 3 . *F ACTOR+ 1 . ) 

IF (ZZ.EQ.5) F F= 1 «, 

8 7 CX = FF/’(CF*CF*DF*DF) 

Z = C X*2 0. 0 / ( PN4QN) 

BT=3.142*DI/QQ-B0 
ZA=BT*BT/ ( 16.0*TS4GC ) 

ZB=0.35*BT/TS 
ZC=HO/ BO 
ZD = H X*. 333/BS 
ZE=HY/BS 

IF ( ZZ-2 ) 88,89,90 

88 PC = ?.4( ZE+ZD+ZA+ZB) 

GO TO 94 

89 PC=Z*< ZC+(2.04HT/(B0+BS) ) + < HW/ BS ) + ZD+Z A + 7.B ) 

GO TO 94 

90 IF (ZZ-4) 91,92,93 

91 PC=Z4(ZC+(2.0*HT/(B0+B1) )+ ( 2 ,0*HW/ ( B1+B2 ) ) + ( HX/ ( 3. *B2 ) I+ZA+ZB) 

GO TO 94 

92 PC=Z*< ZC+0.62) 

GO TO 94 

93 PC=Z*( ZE+Z0+(0.5*GC/TS)+ (0.2 5*TS/GC 1+0.6) 

94 EK=EL/(10.0**(0.103*YY*TS+0.402) > 

IF (DI-8.0) 95,95,96 

95 EK=SQRT( EK ) 

96 ZF=.612*AL0G( 10. 04CS ) 

EW=6.284EK4ZF4< TP**(0.62-< . 2284AL0G( ZF ) ) ) )/( (CL-HV*BV)*DF4DF ) 

C 

C STATOR WINDING LEAKAGE AND ARMATURE REACTION REACTANCES 

C 

IF (CQ) 97,97,98 

97 AB=3.1416*PE 

CQ=(4.*PE+1. )/5. -SIN ( AB)/3. 1416 

98 XL=XR4(PC+EW) 

XD= ( 0. 454EC*PI*CM*0F4100.) /( FH*PX ) 

XQ=( (CQ) / ( CM*C1 ) )*XD 
C 

C POLE LEAKAGE PERMEANCES 

C 

AA=TP*(1.0-PE) 

> X X= PHH+GC 
P1=6.384PHL4XXX/ AA 

P2=4.06*XXX«AL0G( 1 . 0+ 3 . 1 41 59*PHW/ ( 2.0+AA) ) 

P3 = 0 . 

P4=0 . 

DO 99 1=1,4 

>XX= (DR-2.0*PHH)-FL0AT(2*I-1 )*PBH/4.0 
XXX=3. 141594XXX/PX-PBW 

P3=P3+0.44PBL4PBH4(8 .C-FLOAT (241-1 ) >/XXX 

99 P4=P4+0.2544P6H*( 8.0-FLOAT (2*1-1 ) ) 4AL0G ( 1 . 0+3 . 14 1 594PBW/ ( 2.04XXX ) ) 
PEFF=2 .0* ( P1+P2 ) +P3+P4 

C 

C FIELD LEAKAGE REACTANCE, SELF INDUCTANCE AND TIME CONSTANT 

C 

STATET=QQ4SC*0F*CF/( 2.*PN*C) 

XF=XD*(1.0-(Cl/CM)/( 2.04CP+(4.04PEFF/CL )/( 3. 141 594 AG) ) ) 

SI= ( PT442 )4PX4PBL4( ( CP4AG43. 14159 ) /2.0+PEFF/CL )4l.0E-08 
TC=SI/FK1 
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c 

C SYNCHRONOUS AND TRANSIENT REACTANCES CALCULATIONS 

C 

XA=XL+XD 

XB=XL+XQ 

XU=XL+XF 

C 

C DAMPER BAR AND SUBTRANSIENT REACTANCES, ETC. 

C 

IF ( IBN.EO.O) GO TO ICO 
A A= . 62 

IF (DD.EQ.O.) AA=H/ ( 3. *B 1 

BD=(HD/W0+AA+0.5)*6.38 

BE=6.38*( PHW-IBN-L.O ) *T6 ) / ( 3 . 0*GE ) 

PDD=(BD+BE)*(P£FF/CL )*COS( ( 8N-1 . 0 ) *TB* I . 5708/T P ) / ( BD+BE + ( PEFF/CL ) ) 

XDD=XR*PDD 

XX=XL+XDO 

PDQ=(HD/WO+AA+0. 5+GC/TP 1*20. *TB/TP 

XOG=XR*PDQ 

XY=XL+XDQ 

GO TO 101 

100 XDD=0. 

XDQ=0. 

XX = XU 
XY^XB 

101 XN=( XX-t-XY 1/2.0 
TA=XN/(628.3*F*RG1) 

TTC=XU*TC/XA 

C 

C WEIGHT CALCULATIONS 

C 

IF (ZZ-3) 102,103,102 

102 Wl = ( (DU+Dl )*(DU-DI 1*3. 1 A 16 1/ A. 

IF (ZZ.NE.A) WI=WI-QQ*(BS*HS-( < HO+O . 5*HT ) * < BS-BO ) ) ) 

IF (ZZ.EQ.A) WI=WI— QQ*(BS*BS*3.1AL6/A.+H0*80) 

GO TO 10A 

103 WI=(0U-HC)*3.1A16*HC 
WI=WI+HS*( (DI+2.AHS) *3.1A16-GQ«B3) 

WI=WI+QQ*< (HO+O. 5*HT ) *( BS-BO ) ) 

10A WI =wl*0. 283*SS 
C 

RC= ( «321*PT*FE*AS)*PX 
C 

WC=. 321*SC*QQ*AC*HM 
C 

IF (WROTOR.NE.O. ) GO TO 105 
WRCORE=. 283*3. 1A16*( 0SH**2-D ISH**2 )/ A.*PBL 
THETA=2.*3.1A16*PE/PX 
ATIP=DR**2*<THETA-SIN(THETA) 1/8. , 

APHEAD = PHW*< PHH~DR*( 1 . O-COS I THET A/ 2 . 1 ) ) 

AP80DY=PBH*P8W 

BETA=ARSIN(PBW/DSH) 

A8ASE=DSH*PBW*( 1 .O-COS ( BETA ) )-DSH**2* ( 2 . 0*BETA-S I N ( 2 . 0*BETA 1 1/0. 
WP0LE=0.283*(PHL*( AT I P + APHEAO ) +PBL*< APBODY+ ABASE ) ) 
WROTOR=WRCORE+PX*WPOLE 
105 WT0TAL=WC+W1+RC+WR0T0R 
RETURN 
END 
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JsTFK - 1.0 


= STFK 1 



STFK = 1.0- 



12. 5xlQ"4 ‘ 
LT 


lZY • IPX * IPN 
I DM" 0 

46 1* ~~ 

1 1 DM = I DM + IZY I 


\IQQ - 1 DM/ 

< 0 \ / >0 

(Fractional slot winding) \ / 


Reduce IQQ/IZY 
to lowest terms 
(■ IIQQ/IZY) 


Calculate DF 
for fractional 
slot winding 

48-49 1 
Check for 
errors 


= 0 

(Integral 
slot winding) 


Calculate DF 
for integral 
slot winding 


= 0 

vLTRI /to 


Is armature / 
\conductoy , 
\round A— 


AC = 0.785*DW*DW*SN1 


DT = AMINKDW, DWI) 
DG = AMAX1 ( DW, DWI) 


?999 \ RK1 / 
<0. 9999 1 


Calculate "D" (area 
correction for 
corner radii) 


\ ZY / ZY> 0 

zy^oY 


AC = (DT’DG - D)°SN1 


, LTR1 / <0.045 


, LTR1 / <0.094 
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Calculate end turn length, stator 
winding length, stator resistance, 
armature current density 


yls field con-/ 

\ ductor/ No 
\round/ — - — 


AS = 0, 785 # RD*RD 


Calculate field 
resistance 


DT - AMIN1 (RT, RD) 
DG = AMAX1 (RT, RD) 
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SUBROUTINE OUTPUT 
C 

COMMON A,AA,AB,AC s ACORE»Al , ALPHAE » AL PHAR , AL PHA S , A POLE » APOLHD , ARC OR 
1E,AS,AT00TH,B,B1»B2, B3,BCL,BCOIL,BG,BK»BN,BO»BPL,BRC,BS»BTL,BV,C,C 
21,CC,CCR,CE,CF,CL,CM,CQ,CP»Dl,DCORE,DD,DF,DI ,OISH f DR , DRCORE , DSH, DU 
3 »DW,DW1,EC ,EE,EL,ENDEXT, EP ,EW, F, FCL, FE,FFL, FGL1, FGML, FGXL, FH,FKL ,F 
4PL,FQ,FRC,FS»FTL,G,GA,GC»GE,GP,H,HC,HC0IL,HD,HM,H0,HS,HT t HV,HW,HX, 
5HY, IBN, IPN, l PX, IQQ, I WF,I ZZ,KSAT,LTR,LTR1,LTS,P1,P2,P3, P4,PBA,PBH, P 
6BL,PBW,PC,PDD,PDQ,PE,PEFF,PF,PFC,PHH,PHL,PHW,PI,PN,PPL,PPRL,PT,PTC 
7H,PX,QN,GQ,RC,RD»RE,RF,RG1 ,RK,RK1,RPM,RR,RS,RT,RY,S,SB,SC,SD,SF,SH 
8, SI , SIGMA ,SK, SN.SNl, S S , STATE T , T 1 , T 1 1 t T2 , T2 2 , T 3 , T 33, T A, TB , TC , TF , TG , 
9TS, TST,TT,TTC,VA,VMI N,VR,WC,WI , WL , WO, WROTOR , WTOT AL , X A, XB , XD , XDD, XD 
$Q,XF,XL,XN,XQ,XR,XU,XX,XY,YY,ZG,ZZ 
C 

DIMENSION G ( 5 ) , PFC ( 5 ) , 41(60) 

C 

INTEGER ZZ 
C 

REAL LTS, LTR , LTRL 


WRITE (6,1) VA,E£»EP»PI» PF» IPN,F, IPX,RPM 

FORMAT (1HL,18H ALTERNATOR R AT I NG/ / LOX , 1 5H ALTERNATOR KVA, F16. 1/10 
1 X , 1 8H LINE-LINE VOLT AGE, F12.0/ LOX, 19H LINE-NEUT. VOL T AGE , F 1 1 . 0/1 OX 
2 , 1 AH PHASE CURRENT, F18.2/10X,13H POWER F ACTOR, F 19 . 2 / 10X , 7H PHASES, 
3I22/10X»10H FREQUENCY, F20.0/10X.6H POLES , I 2 3/ 10X , 4H RPM,F27.1) 

IF (IZZ-?) 3,5,2 

IF (IZZ-4) 7,9,11 

WRITE (6,4) BS,HX,HY,HS, IQQ, TS,TT 

FORMAT ( 1 HL , 1 3H STATOR SL0TS//5X10H TYPE-OPEN/ 54X , 9H *,12X6 

1H* /62X1H*, 12X1H*/55X2HHY,5X1H*, 12X1H*/10X3H BS.F26.3, 1 X6H INCH 

2ES, 16X, 1H*,12XIH*/10X3H HX , F 26 . 3 , 15X , 9H *, 2X8H********, 2X1H 

3*/ 1 0 X3H HY,F26.3,23X1H*, 2X1H*,6X1H*, 2X1H*/10X3H HS , F 26 . 3 , 2 3X 1H*, 2X 
41H*,6X1H*,2X1H*/62X» 1 H * , 2X8H********, 2X 1 H*2X 2HHS/55X2HHX , 5X , 1 H*, 12 
5X1H*/10X13H NO. OF S LOTS 1 16 , 23X, 1 H*, 2X8H********, 2X 1 H*/62X 1 H*, 2X 1H 
6*,6X1H*,2X1H*/10X11H SLOT PITCH,F18.3, 1X6HINCHES , 16X 1H* , 2X1 H*, 6X 1H 


7*,2X1H*/54X9H * , 2 X8 H* ******* , 2X1H*/10X11H SLOT PITCH, 41X1H* 

8, 12X1H*/10X15H AT 1/3 DIST.,F14.3,1X6HINCHES,16X 1 9 H** *********** 


9* /62X1H1, 12X1H1/62X14H1 BS 1/62X1 HI, 12X1H1 ) 

GO TO 13 

WRITE (6,6) BO,BS,HO,HX,HT,HW,HS, IQQ,TS,TT 

FORMAT (1HL.13H STATOR SL0TS//5X22H TYPE-PART I ALLY CL0SED/67X4H-B0 

1-/5 7X1 OH *»4X10H* /5 8X2HH0 ,6X1H*,4X1H*/57X10H 

2 *, 4X1H*/ 10X3H BO, F26.3, 1X6HINCHES, 19X1H*,6X1H*/10X3H BS,F26. 3, 

3 19X2HHT ,4X1H*,8X1H*/10X3H HO, F26. 3, 24X1H*, 10X1H*/10X3H HX,F26.3,18 

4 X6H *, 12X1H*/10X3F HT, F26. 3, 23X1H*, I2X 1H*/ 10X3H HW, F 26 . 3 , 1 9X2H 

5HW , 2 XI H* , 1?X1H*/10X3H HS , F26 . 3 , 1 8X6H *2X8H********, 2X1H*, 2X2HH 

6S/62XIH*, 2 XI H* ,6X1H*,2X1H*/10X13H NO. OF SLOTS 1 16 , 2 3X 1H* , 2X 1 H* , 6 X 1 
7H*, 2X1 H*/62X1H*,2X8H* *******, 2 X1H*/ 10X1 1H SLOT P I TCH , F 1 8 . 3 , 1 X6HI NC 
8HES, 12 X2HHX,2X1H*, 12X1 H*/62X1H*,2X'8H********,2X1H*/ 10X1 1H SLOT PIT 
9CH,41X1H*,2XIH*,6X1H*,2X1H*/ 10X15H AT 1/3 DIST.,F14.3, 1X6HINCHES 

$ ,16X1H*2X1H*6X1H*2X1H*/5 7X6H *, 2X8H********, 2X1H*/62X1H*, 12X 1H 

$*/ 62 X19H* ************* /62X1H1, 1 2X 1H 1/ 62 X 1 4H1 BS 1/62X 

$ 1H1 , 12X1H1 ) 

GO TO 13 

7 WRITE (6,8) BO,B1,B2,B3,BS,HO,HX,HT,HW,HS,IQQ,TS,TT 
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8 FORMAT ( 1 H L » 1 3 H STATOR S LOTS //5X25HTYPE-CONST ANT TOOTH WIDTH/ 6LX1H 

1 1 , 14X1H1 / 61X16H1 B 1 1/10X3H BO , F 26 . 3 , 1X6H I NC HE S, 15X1H1 , 1 

24X1H1/10X3H B1,F26.3 »22X1H1, 5X4H-B0-»5X1H1/10X3H B2 , F26. 3 » 1 IX 1 7H— 

3 — — — 1- — -*, 4X17H* 1 — — / 1 OX 3H B3.F26.3, 22X1H1,4X1H* 

4,4X1H*,4X1H1 , 8 X2HHO/ ICX15H BS = (B2+B3)/2,F14.3,22X1H1,4X1H*,4X17H 

5*-- — 1 / 10 X3H HO, F 26. 3, 22X 1H1 , 2X 1H* , 8 X 1H* , 2X 1H1 , 8X2HHT / 1 

60X3H HX,F26.3,22X1H*,14X,12H*— /10X3H HT, F26. 3, 12X2HHS, 7X 

7 1H* , 16X1H*,7X2HHW/10X3H HW , F26. 3 , 20X 1H*, 3X 1 2H************ , 3X 1 OH* — 

8 /10X3H HS,F26.3,19X2H*1,3X1H*, 10XIH*, 3X2H1*/57X 1H* , 1 X1H1 , 3X 

9 1H* , 10X1H*,3X1H1, 1X1 H* ,4X2 HHX/ 10X 1 3H NO. OF SLOTS , 1 16, 17X1H* , 2X1 HI 

$,3X12H*** *********, 3X1 HI ,2X1H*»6H / 55X 1 H* , 3X1H1 , 1 8X 1H 1 , 3X1H*/ 

S10X11H SLOT PITCH»F18.3,1X6HINCHES» 4X34H ********* ****** ******* 

$********/5AXlHl,4XlH 1,18X1 HI, 4X1H1/10X11H SLOT PITCH, 33X1H1»4X20H1 

$ B 2 1 , 4X 1H1 / 10XI5H AT 1/3 DIST. ,F14.3, 1X6H INCHES, 8 

tXIHl ,4X1H1,18X1H1 ,4X1H1/54X30H1 B3 1/54X1H 

$1,28 XI HI) 

GO TO 13 

9 WRITE (6,10) BO,HO,BS,HS,IQQ,TS,TT 

10 FORMAT (1HL.13H STATOR S LOTS //5X , 1 1 H TYP E-ROUND/ / 10X , 1 3H SLOT OPEN 
1 ING,F16.3, IX 6H INCHES / 10X,19H SLOT OPENING DEPTH, F 10 . 3/ 10X , 1 4H SLOT 
2 DIAMETER, F15.3/10XUH SLOT DEPTH, F18. 3//10X, 13H NO. OF SLOTS, 116/ 
3/10X.11H SLOT PITCH, F18.3, 1X6HINCHES//10X, 11H SLOT P ITCH/10X, 15H 

4 AT 1/3 DIST. ,F14.3, 1X6HINCHES) 

GO TO 13 

11 WRITE (6,12) BS,HX,HY,HS,IQQ,TS,TT 

12 FORMAT (1HL.13H STATOR SL0TS//5X25H TYPE-OPEN (1 COND. /SLOT ) /57X , 6 

1H *1 2 X 6 H* /62X, 1H*,12X1H*/58X5HHY * , 12X lH*/62 X1H* , 1 2X 1 H*/ 

2 10X, 3H BS,F26.3, 1X6H INCHES, 11X.6H * , 2X8H******** , 2X1 H*/ 1 OX, 3H 

3HX,F26.3,23X,1H*,2X1H*,6X1H*,2X1H*/10X, 3H HY, F26. 3, 23X, 1H*, 2X1H* ,6 
4X1H*,2X1H*/10X,3H HS , F26 .3 , 2 3X , 1H* , 2X 1H*, 6X 1H* , 2X 1H*, 2 X2HHS/58X2 HH 
5X,2X1H*,2X1H*,6X1H*,2X1H*/10X,13H NO. OF SLOTS , 1 16, 23X 1H* , 2X 1H*, 6X 
61H*,2X1H*/62X1H*,2X1H*,6X1H*,2X1H*/10X,11H SLOT PITCH, FI 8 . 3, 1X6HIN 

7CHES,16X1H*,2X1H*,6X1H*,2X1H*/57X6H *, 2X8H********, 2X1H*/10X 11 

8 H SLOT PITCH, 41X1H*, 1 2X1 H*/ 1 OX 1 5H AT 1/3 D I ST . , F 1 4 . 3 , 1 X6H I NCHE S , 

9 16X1 9H*** *********** /62X1H1, 12X1H1/62X14H1- BS 1/62X1H 

$1 , 12X1H1 ) 

13 CONTINUE 
C 

WRITE (6,14) GC , GP , GE 

14 FORMAT ( 1HL , 8 H AIR GAP//10X,16H MINIMUM AIR GAP, F 1 7 . 3, 1X6H INCHES / 1 
10X.16H MAXIMUM AIR GAP ,F17.3/10X, 18H EFFECTIVE AIR GAP.F15.3//) 

IF ( IBN.EQ.O) GO TO 16 
WRITE (6,15) CC.CCR 

15 FORMAT tlH , 10X, 18HCARTER COEFFICI ENT/17X, 6HSTAT0R, F20. 3/17X5HR0T0 
1R, F21. 3) 

GO TO 18 

16 WRITE (6,17) CC 

17 FORMAT (1H , 10X, 18HC ARTER COEFF I C I ENT , F 1 4. 3 ) 

18 CONTINUE 

IF (RF.LT..5) WRITE (6,19) 

IF (RF.GE..5) WRITE (6,20) 

19 FORMAT (IH1.45H ARMATURE WINDING ( Y-CONNECTED, RANDOM WOUND)//) 

20 FORMAT (1H1.43H ARMATURE WINDING (Y-CONNECTED, FORM WOUND)//) 

IF (DW1.EQ.0. ) WRITE (6,21) DW 

21 FORMAT (1H ,9X,16H STRAND D I AMETER, F 32. 4, 1 X6H I NCHE S ) 

IF (DW1.GT.0.) WRITE (6,22) DW,DW1 

22 FORMAT ( 1H ,9X,18H STRAND D I MENS I ONS , F30 . 4 , 2H X, IX, F 6 . 4, IX, 6 HINCHE 
IS) 
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WRITE (6,23) SD,SN,SN1,AC,S,CE,HM,EL,ENDEXT,SK,RS,TST,RG1,$TATET,E 
1C,PTCH,YY,QN 

23 FORMAT (1H » 10X36HDI STANCE BTWN CL OF STRANDS ( R AD I AL ) , F 1 1 . 4// 10X , 
123H STRANDS/CONDUCTOR IN RADIAL DIR. , F 1 1 . 0/ 1 OX , 24H TOTAL STRANDS/C 
2CNDUCTOR,F20.0/10X,15H CONDUCTOR AREA, F 3 3. 4, IX 6HSQ- I N . / L OX , 29H CUR 
3RENT DENSITY AT FULL L OA D, F 17 . 2 , 3X 1 OH AMP/ SQ- I N . / / 1 OX , 27H COIL EXTE 
4NSI0N BEYOND CORE, F20. 3, 2X6H INCHES/ 10X, 2*H MEAN LENGTH OF 1/2 TURN 
5,F23„3/10X,16H END TURN LE NGTH , F 31 , 3/10X 19H END TURN EXTENSI 0N,F28 
6.3/10X.17H STATOR SLOT S KEW , F30 . 3// 1 OX, 25H RESISTIVITY AT 20 DEG. 
7C.F23.4, 1X16HMICR0 OHM I NCHE S/ 10X , 2 1H STATOR RESISTANCE AT,F6.0,7H 
8 DEG. C*F14. 4, 1X4HOHMS//IOX, 3 OH NO. OF EFFECTIVE SERIES TURNS, F16. 
92/10X.27H TOTAL EFFECTIVE CONDUCTORS , F 1 9 . 2/ 1 OX , 1 4H WINDING PITCH, F 
$34. 4/10X, 14H SLOTS S P ANN ED, F 30 . 0/ 10X , 25H SLOTS PER POLE PER PHASE, 
SF21.2) 

WRITE 16,24) SC,C,PBA,FS»DF,CF 

24 FORMAT (1H ,9X,16H C ONDUCTOR S/ SLOT , F28. 0/ 10X , 25H NO. OF PARALLEL C 
1IRCUITS,F19.0/10X»17H PHASE BELT ANGLE , F 27 . 0 , 5X7HDEGREES // 10X , 12H 
2SKEW FACTOR, F35. 3/ 10 X , 20H DISTRIBUTION FACTOR, F27 . 3/ 10X, 13H PITCH 
3F AC TOR , F34 . 3 ) 

IF (RT.EQ.O. ) WRITE (6,25) RD 
IF (RT.GT.O. ) WRITE (6,26) RD, RT 

25 FORMAT (1HL,14H FIELD WINDING//10X19H CONDUCTOR DI AMETER , F29 . 4, 1 X6 
1HINCHES/) 

26 FORMAT (1HL,14H FIELD W I ND I NG/ / 10X , 2 1H CONDUCTOR D I MENS I ON S , F27. 4, 
12H X, 1 XF6. 4, 1X6H INCHES/ ) 

WRITE (6,27) AS,PT,FE,RR,TF»FK1 

27 FORMAT (1H ,9X15H CONDUCTOR ARE A, F 33 . 4, 1 X6HSQ- IN . / / 10X , 2 4H NO. OF 

1 TURNS (PER POLE) , F20 „ 0/ 1 0X20H MEAN LENGTH OF TURN, F27. 3, 2X6HINCHES 
2//10X25H RESISTIVITY AT 20 DEG. C, F23.4, 1X16HM ICRO OHM INCHES/10X2 
30H FIELD RESISTANCE AT,F5.0,7H DEG. C , FI 6. 4, 1X4H0HMS/ 1 7X 1 8H (COILS 
4 IN SERIES)) 

WRITE (6,28) HCOI L,BCO IL 

28 FORMAT (10X12H COIL HE IGHT , F 35 . 3/ 10X 1 1H COIL WIDTH,F36.3) 

WRITE (6,29) Cl ,CP,CM,CQ,D1 

29 FORMAT (1HL.10H CONS T ANTS/ / 1 OX , 35H Cl, FUNDAMENTAL /MAX . OF FIELD F 
1LUX»F8.3/10X» 18H CP, POLE CONSTANT , F25 . 3/ 10X , 2 7H CM, DEMAGNETI ZAT I 
2CN FACTOR, F16. 3/10X, 31.H CQ, CROSS MAGNETIZATION F ACTOR , F 12 . 3/ 10X , 2 
36H Dl, POLE FACE LOSS FACTOR , F17 . 3 ) 

WRITE (6,30) DI ,DU»CL»SS»HC,SF,HV,BV»BK,WL 

30 FORMAT <1H1,7H STATOR// 1 OX , 2 3H STATOR INSIDE D I AMETER , F2 1 . 2 , 1 X6H I N 
1CHES/10X, 24H STATOR OUTSIDE D I AMETER , F20 . 2/ 1 OX , 20H OVERALL CORE LE 
2NGTH,F24.2/10X,22H EFFECTIVE CORE LENGTH, F22 . 2/ 10X , 17H DEPTH BELOW 

3 SL0T,F27«2//10X,16H STACKING FACTOR, F28 . 2// 10X , 21H NO. OF COOLING 

4 DUCTS, F21.0/10X,15H WIDTH OF DUCTS, F29. 2, 1X6H INCHES// 10X, 13H CORE 

5 LOSS AT , F6„ 1 » 1 7H KI LOL I NE S/ SQ . I N . , F7 . 1 , 2X9HWATTS/LB. ) 

IF (LTS.NE.O. ) WRITE (6,31 ) LTS 

31 FORMAT ( 10X, 21H LAMINATION THICKNESS, F24. 3, 4H IN.) 

WRITE (6,32) PBW , PBL , PBH ,RK 

32 FORMAT (1HL,6H R0T0R//10X, 16H POLE BODY WI DTH, F24. 3, 7H I NCHES/20X, 
113H AXIAL LENGTH, F 17 „ 3/ 20X , 7H HE IG'HT, F23. 3/20X, 16H STACKING FACTOR 
2 , F 14 . 3 ) 

IF (LTR.NE.O.) WRITE (6,33) LTR 

33 FORMAT (IH , 19X,21H LAMINATION THICKNESS, F9. 3, 7H INCHES) 

WRITE (6,34) PHW»PHL» PHH,RK1 

34 FORMAT (1HK,9X,16H POLE HEAD WIDTH, F24. 3, 7H INCHES/20X, 1 3H AXIAL L 
1ENGTH,F17.3/20X,7H HE IGHT, F23. 3/20X, 16H STACKING FACTOR, F 14. 3 ) 

IF ( LTR1.NE.0. > WRITE (6,33) LTR1 
WRITE (6,35) PE, DR, VR, SIGMA 
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FORMAT ( 1HK,9X, 13H POLE EMBRACE , F27. 3/10X, 15H ROTOR D I AME TER , F25 . 3 
1/lOX ,17H PERIPHERAL S PEE D, F20. 0 , 3X , 10H FEET /M I N . / 10X f 2 3H SPEC, TAN 
2GENTIAL FORCE , F 1 7 . 3 » 1 1 H LBS/SO. IN.) 

WRITE (6,36) DSH.DISH 

36 FORMAT ( 1HK,9X,23H DIAMETER OF ROTOR CORE, F18.3, 7H INCHES/lOX, 3AH 
1INSIDE DIAMETER (OF HCLLOW SHAFT ) » F7. 3 ) 

IF (IBN.EQ.O) WRITE (6,37) 

37 FORMAT (1HL,19H DAMPER BARS (NONE)) 

IF (DD.EQ.O..AND.IBN.NE.O) WRITE (6,38) H,B 

38 FORMAT (1HL,26H DAMPER BARS ( RECTANGULAR )// 10X , 22H DAMPER BAR DIME 
1NSI0NS,F17.3,2H X, 1XF 5 . 3 , 1 X6HI NCHES ) 

IF (DD.NE.O..AND.IBN.NE.O) WRITE (6,39) DD 

39 FORMAT (1HL,20H DAMPER BARS ( ROUND )// 10X , 20H DAMPER BAR DIAMETER, F 
119.3,1X6HINCHES) 

IF (IBN.NE.O) WRITE (6, AO) WO, HD, SB, TB, I BN, RE, T3 

AO FORMAT ( IH ,9X,I9H SLOT OPENING WIDTH, F20. 3/10X, 20H SLOT OPENING H 
1EIGHT,FI9.3/10X,18H DAMPER BAR LENGTH, F21 . 3/ 10X , 17H DAMPER BAR PIT 
2CH,F22.3//10X»2AH NO, OF DAMPER B ARS/ POL E , 1 1 2// IOX , 2 5H RESISTIVITY 
3 AT 20 DEG. C,FIA.3,17H MICRO-OHM INCHES//IOX, I8H TEMPERATURE (HOT 
A),F20.2,2X5HDEG.C) 

WRITE ( 6 , A1 ) WC,RC,WI , WROTOR , WTOT AL 

A 1 FORMAT (IH1.8H WEIGHTS// 10XI3H STATOR COND. , F17. 3, IX6HP0UNDS/10X12 

1H FIELD COND. ,F18.3/ 10X12H STATOR I RON, F 18 . 3/ IOX, 6H ROTOR, F2A. 3//1 
2 OX , 6H TOT AL/ 1 1 XI 8H ( ELECTROMAGNET I C >, Fll . 3 ) 

WRITE <6,A2) PC, EW 

A2 FORMAT (1HL,I1H PERMEANCES// IOX, 16H WINDING LEAKAGE/10X, 26H (PER I 
INCH OF CORE LENGTH )/ 1 2X, 12H STATOR SLOT , F26 . 3 , 1 8H LINES/AMPERE TUR 
2N,/13X,10HSTATOR END,F27.3) 

IF (IBN.NE.O) WRITE (6, A3) PDO,PDQ 

A3 FORMAT (1HK,9X25H DAMPER-DIRECT (PER INCH) , F15.3/10X, I8H DAMPER-QU 
IADRATURE,F22.3) 

WRITE ( 6 , A A ) P1,P2,P3,PA,PEFF 

AA FORMAT (IHL,9X19H PI POLE HEAD I NNER, F2 1 . 3/ IOX 17H P2 POLE HEAD END 
1 ,F23.3/10X19H P3 POLE BODY I NNER , F2 1 . 3/ IOX 17H PA POLE BODY END,F23 
2.3/10X23H EFFECTIVE POLE L EAKAGE, F 17. 3 ) 

WRITE ( 6 , A5 ) A,XR,XL,XD,XQ»XA,XB,XF,XU 

A 5 FORMAT (1HL,11H REACTANCES//10X23H AMPERE CONDUCTORS/ I NCH, F20. 3/ 10 

1X17H REACTANCE F ACTO R , F26. 3/ /10X23H STATOR WINDING LEAKAGE ,F20. 3, 1 
2X,7HPERCENT/10X23H ARM. REACTION (DIRECT), F20.3/10X22H ARM. REACTI 
3CN (QUAD. ) ,F21.3/10X21H SYNCHRONOUS (DIRECT), F22.3/10X20H SYNCHRON 
AOUS (QUAD. ),F23.3/10X1AH FIELD LEAKAGE , F 29. 3/ IOX 10H TR ANS I ENT, F33 . 
53) 

IF (IBN.NE.O) WRITE (6,A6) XDD, XDQ 

A6 FORMAT (1H ,9X2AH DAMPER LEAKAGE ( D I RECT ) , F19 . 3/ IOX , 23H DAMPER LEA 
1KAGE (QUAD. ) ,F20.3) 

WRITE { 6 , A7 ) XX,XY,XN,SI ,TF,TC*TTC,TST,TA 

A 7 FORMAT (1H ,9X22H SUBTRANSIENT ( D I RECT ) , F2 1 . 3/ 10X2 1H SUBTRANSIENT 

1 (QUAD. ),F22.3/10X18H NEGATIVE SEQUENCE, F25. 3// 10X22H FIELD SELF IN 
2DUCTANCE,F21.3,1X7HHENRIES//10X27H OPEN CIRCUIT TIME C0NSTANT/13X, 
3 16H (FIELD ONLY, AT,FA.0,8H DEG. C ) , F 1 A . 5, 1X7HSEC0NDS/ 10X2AH TRANS 
A I ENT TIME CONSTANT, F21 .5/10X23H ARMATURE TIME C0NSTANT/13X12H (WIN 
5DING AT,FA.Q,8H DEG. C),F18.5) 

RETURN 

END 
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SUBROUTINE MAGNET 


COMMON a,aa,ab,ac,acore, ai » alphae, alphar, alphas, apole.apolhd, arc or 
1E,AS,AT00TH,B,B1,B2, B3,BCL,BCOIL,BG,BK,BN, BO , BPL , BRC , B S, BTL , BV ,C , C 
21,CC,CCR,CE,CF,CL,CM,CQ, CP , D I , DCORE, OD , DF, 0 1 , D I SH, DR , DRCORE , D SH, DU 
3,DW,DW1,EC,EE,EL,ENDEXT,EP,EW,F,FCL,FE,FFL,FGL1,FGML,FGXL,FH,FK1,F 
APL,FQ,FRC,FS,FTL,G,GA,GC,GE,GP,H,HC,HCOIL,HD,HM,HO,HS,HT,HV,HW,HX, 
5HV,IBN,IPN,IPX,IQQ,IWF, I ZZ , KSAT , LTR, LTR 1 , LTS , P l , P2 , P3, PA, PB A , PRH , P 
6BL,PBW,PC,PDD,PD0,PE , P EFF, PF , PFC , PHH, PHL , PHW , P I , PN, PPL , PPRL , P T , P TC 
7H,PX,QN,QQ»RC,RD»RE»RF»RG1»RK»RK1,RPM»RR,RS»RT,RY,S,SB»SC»SD»SF,SH 
8, SI ,SIGMA,SK,SN,SN1,SS,STATET,T1,T11,T2,T22,T3,T33,TA,TB,TC,TF,TG, 
9TS, TST.TT.TTC, VA.VMI N, VR.WC, WI , WL , WO , WROTOR , WTOTAL , X A, XB,XD,XDD,XD 
$Q,XF,XL,XN,XG,XR,XU, XX, XY, YY,ZG,ZZ 
C 

DIMENSION G { 5 ) , PFC( 5 ) , AI(60) 

C 

INTEGER ZZ 
C 

REAL LTS,LTR,LTR1 
C 

BPL=0. 

BTL=0. 

BCL=0. 

BRC =0. 

FFL=0. 

FPL=0. 

FCL =0. 

FTL=0. 

FRC=0. 

PPRL=0. 

FAGT=FGH+FGXL 

FLUX DENSITY AND AMPERE-TURNS FUR TEETH 


BTL=PPL/ATOOTH 
X=BTL 
NA= 1 
K= 1 

KSAT=3 
GO TO 6 
FTL=AT*HS 

FLUX DENSITY AND AMPERE-TURNS FOR CORE 

BCL=PPL/{ 2.0#AC0RE) 

X=BCL 
NA = 1 
K SA T=4 
K=2 

GO TO 6 
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FCL= AT*DCORE 


POLE HEAD LEAKAGE FLUX 

FPRL = ( FAGT + FTL + FCU* ( 2.0*( P1+P2 S> *0.001 

FLUX DENSITY AND AMPERE-TURNS FOR ROTOR POLE HEAD 

BPL= (PPL+PPRL) /APOLHD 
X=BPL 
NA = 3 1 
KS AT=1 
K = 3 

GO TO 6 
FPL = AT *PHH 

FLUX DENSITY AND AMPERE-TURNS FOR ROTOR POLE BODY 

XXX=FAGT+FTL+FCL+FPL 
PPRL=PPRL+XXX*(P3+P4)*.001 
X=( PPL+PPRLI/APOLE 
BPL = AMAX1 ( BPL , X ) 

KSA T = 1 
K=5 
NA = 3 1 
GO TO 6 

FPL=FPL+AT*PBH 

FLUX DENSITY AND AMPERE-TURNS FOR ROTOR CORE 

BRC= I PPL+ PPRL)/(2.0*ARCGRE) 

X = BRC 
NA = 3 1 
KSA T=2 
K = 4 

GO TO 6 
FRC=AT*0RC0RE 

TOTAL AMPERE-TURNS 

FFL=FAGT+FTL+FCL+FPL+FRC 
KS A T=1 0 
RETURN 

INTERPOLATION PROCEDURE FOR MATERIAL CURVES 

IF ( AI (NA).LT.X) GO TC 10 
NA=NA+3 

7 IF (AI(NA)-X) 8,9,9 

8 NA=N A+2 
GO TO 7 

9 >X=< A I (NA)-AI (NA-2I) /< ALOGI AKNA+L )/( A I (NA-1 )+0.0001) ) ) 
Y = A I (NA) — XX* ALOGI AKNA + l I) 

AT= E XP ( I X— Y ) /XX ) 

GO TO 11,2,3,5,4) ,K 

10 RETURN 
END 
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SUBROUTINE WINDGE ( G , RPM ,DR , CL , KW I ND, GC , N, WFTTL ) 

DIMENSION G( 5 ) , PRESS ( 5 ) * FLUID<6) 

REAL M 

NAMELIST /WIND/ VIS, PRESS, PRSRAT,M, TEMP, RO 

IF (KWIND.NE.i) GO TO 3 

R0=0 

N*1 

DO I 1=1,5 

1 PRESS(I)=0. 

READ ( 5,2) FLUID 

2 FORMAT (6A6) 

READ (5, WIND) 

PRSR=PRSRAT 

IF (RO.EQ.O.) R0=(PRSR*M*5.l8E-02)/(TEMP+273.0) 

GO TO 4 

3 IF (PRESS(N).EQ.O.) PRESS! N ) =PRSRAT*G(N ) 

PRSR=PRESS(N> 

RO=( PRSR*M*5.18E-02) / (TEMP+273.0) 

MAIN AIR GAP 

RE 1 = 3. 14l59*(DR/2.0) *GC*RPM*RO/ <30.* 144.4V IS) 

A=AL0G(RE1 ) 

IF (RE1.LT.7000. ) CD= 10. 59346-1 ,78872*A+.08205*A**2 
IF (RE1.GE.7000. ) CD=15. 45399-2. 81449*A+. 135502*A**2 
CD=EXP(CD)*1.0E-03 
IF (RE1. GT. 32000. ) C0=2.25E-03 

WF1 = (6.107E-10)*CD*R0*(RPM**3 )*( ( DR/2.0 ) **4)*CL 
WF T TL= l . 5*WF 1 
IF (KWIND.NE.I) RETURN 

WRITE (6,5) FLUID, VI S,RO,PRSRAT, TEMP, M 

FORMAT ( 1H1 , 22H WINDAGE AT RATED L0AD//10X20H FLUID PROPERTIES OF, 
11X6A6//13X.10H VISCOSITY, 1PE25. 3, 11H LBS/SEC-FT/ 13X8H DENS I TYOPF 24 
2.4, UH LBS/CU.FT./13X.9H PRESSURE, 0PF21.2, 11H LBS/SQ. IN/13X,12H T 
3EMPERATURE,F18.2,7H DEG. C/13X,17H MOLECULAR W E I GHT , F 1 3 . 2 ) 

WRITE (6,6) RE1 ,WFTTL 

FORMAT ( 1HK »9X , 16H REYNOLDS NUMBER, F 17 . 2// 10X1 3H WINDAGE L0SS,F20. 
12,1 X5HWAT T S ) 

RETURN 
END 
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APPENDIX D 


DEFINITION OF FORTRAN VARIABLES 


The following is an alphabetic list of the FORTRAN variables used in the program. 
Each variable is defined, and the units used in the program for each variable are given. 
Variables that are subscripted (dimensional) or that are permissible program input are 
identified by a + or *, respectively. In addition to some of the previous figures in- 
cluded in the report, figures 17 and 18 are used in the definitions of the FORTRAN vari- 
ables. 


A 

A 

AA 

AB 

.ABASE 

AC 

AC ORE 
AG 
*+ AI 

AIRGAP 
-f‘ A KVA 
ALPHA 

* ALPHAE 

* ALPHAR 

* ALPHAS 
AN 

APBODY 

APHEAD 

APOLE 

APOLHD 

ARC ORE 


ampere conductors per inch of stator periphery, A/in. 

dummy variable used in subroutine WINDGE 

dummy variable used for variety of calculations 

dummy variable used for variety of calculations 

o 

area used in pole weight calculations, in. 

2 

armature conductor area, in. 

2 

stator core cross-sectional area used in magnetic calculations, in. 

specific air gap permeance per pole, lines/(A-turn)(in. ) 

coordinates of points on material magnetization curves (see figs. 9 
and 17) 

NAMELIST name 

alternator output at load G, kVA 

angle used in end-turn length calculations (see fig. 18), rad 
temperature coefficient of RE, per °C 
temperature coefficient of RR, per °C 
temperature coefficient of RS, per °C 
power factor angle (see fig. 4), rad 

2 

area used in pole weight calculations, in. 

2 

area used in pole weight calculations, in. 

2 

pole area used in magnetic calculations, in. 

2 

pole head area used in magnetic calculations, in. 

2 

rotor core cross-sectional area used in magnetic calculations, in. 
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AS 

o 

field conductor area, in. 


AT 

ampere-turns per inch at flux density X for magnetic materials (see 
fig. 17), A-turn/in. 


ATIP 

2 

area used in pole weight calculations, in. 


ATOOTH 

o 

tooth area used in magnetic calculations, in. 

* 

B 

rectangular damper bar width, in. 


B1 

stator slot dimension (see fig. 10), in. 


B2 

stator slot dimension (see fig. 10), in. 

* 

B3 

stator slot dimension (see fig. 10), in. 


BARD 

damper bar equivalent diameter used in loss calculations, in. 


BB 

displacement angle (see fig. 4), rad 


BCL 

2 

stator core flux density, kilolines/in. 

+ 

BCLL 

2 

stator core flux density, kilolines/in. 

* 

BCOIL 

field coil width (see fig. 14), in. 


BD 

used in damper bar leakage permeance calculations 


BE 

used in damper bar leakage permeance calculations 


BETA 

angle used in pole weight calculations, rad 


BG 

2 

main air-gap flux density (no-load rated voltage), kilolines/in. 

* 

BK 

p 

flux density at which core loss WL is specified, kilolines/in. 

* 

BN 

number of damper bars per pole 

* 

BO 

stator slot dimension (see fig. 10), in. 


BPL 

2 

pole flux density, kilolines/in. 

+ 

BPLL 

2 

pole flux density, kilolines/in. 


BRC 

2 

rotor core flux density, kilolines/in. 

+ 

BRCORL 

2 

rotor core flux density, kilolines/in. 

* 

BS 

stator slot dimension (see fig. 10), in. 


BT 

dummy variable used in slot and end turn leakage permeance calcula- 
tions 


BTL 

2 

tooth flux density, kilolines/in. 

+ 

BTLL 

2 

tooth flux density, kilolines/in. 
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* BV 

* BX 

* C 

* Cl 
CALPHA 
CC 
CCR 

CD 

+ CDD 

* CE 
CF 
CK 

* CL 

* CM 

+ CMPNT 

* CP 

* CQ 
CS 
CW 

cx 

CY 

CZ 

D 

D 

* D1 
D2 
D3 
D4 
D5 
D6 


width of cooling ducts, in. 

thickness of armature coil at end-turn (see fig. 11), in. 
number of parallel armature circuits per phase 

ratio of fundamental maximum to actual maximum value of field form 
cosine of alpha 

Carter coefficient for stator slots 

Carter coefficient for rotor slots 

skin friction coefficient of gas in alternator air gap 

2 

current density in field winding, A/in. 
armature coil extension (see fig. 11), in. 
pitch factor 

power-factor adjustment factor 

length of stator stack (axial direction), in. 

demagnetizing factor (direct axis) 

array of names used in saturation message printout 

ratio of average to maximum value of field form 

cross-magnetizing factor (quadrature axis) 

armature winding pitch, per unit 

winding constant 

used in armature conductor slot leakage permeance calculations 

2 

used to calculate area correction D, in. 

2 

used to calculate area correction D, in. 
used in distribution factor calculation 

2 

area correction for corner radii in rectangular conductors, in. 

pole face loss factor 

factor used in pole-face loss calculation 

factor used in pole-face loss calculation 

factor used in pole-face loss calculation 

factor used in pole-face loss calculation 

factor used in pole-face loss calculation 
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DA 


* 


* 

* 


+ 


* 


* 

* 

* 

+ 

+ 

+ 

+ 

+ 


* 


+ 


* 


DAMPER 

DB 

DC ORE 

DD 

DF 

DG 

DI 

DISH 

DL 

DR 

DRCORE 

DSH 

DT 

DU 

DW 

DW1 

DX 

DY 

DZ 

E 

El 

EB 

EC 

EDD 

EE 

EF 

EK 

EL 


2 

used to calculate area correction D, in. 

NAMELIST name 

diameter of bender pin for forming armature coils, in. 
stator core length used in magnetic calculations, in. 
damper bar diameter, in. 
distribution factor 

maximum dimension of rectangular conductor, in. 
stator inside diameter, in. 

rotor inside diameter (for rotor with center hole) (see fig. 14), in. 

damper bar loss, W 

rotor diameter at main air gap, in. 

rotor core length used in magnetic calculations, in. 

outside diameter of rotor core (see fig. 14), in. 

minimum dimension of rectangular conductors, in. 

stator outside diameter, in. 

armature conductor strand diameter or width (see fig. 11), in. 

armature conductor strand height (uninsulated) (see fig. 11), in. 

2 

used in rectangular conductor area calculation, in. 

2 

used in rectangular conductor area calculation, in. 

2 

used in rectangular conductor area calculation, in. 
alternator electromagnetic efficiency, percent 
alternator overall efficiency, percent 
eddy factor (bottom) 

number of effective armature conductors 
excitation voltage, per unit 
line-to-line rated alternator voltage, V rms 
field voltage per coil, V dc 

leakage reactance factor for armature winding end-turns 
length of end-turn for armature coil, in. 
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END EXT 


EP 

ET 

EW 

EX 

EZ 

F 

FACTOR 

FAGT 

FCL 

FCLL 

FDMLL 

FE 

FF 

FFL 

FFLL 

FGLl 

FGLL1 

FGML 

FGXL 

FH 

FI 

FIELD 

FIMM 

FK1 

FLUID 

FNQ 

FPL 


horizontal extension of armature coil end-turn beyond stator core (see 
fig. 18), in. 

line-to-neutral rated alternator voltage, V rms 
eddy factor (top) 

specific armature conductor end-turn leakage permeance, 
lines/(A-turn)(in. ) 

eddy loss, W 

eddy factor (average) 

frequency, Hz 

used in armature conductor slot leakage permeance calculations 

FGLl + FGXL, A-turn 

stator core ampere turns, A-turn 

stator core ampere turns, A-turn 

direct-axis demagnetizing ampere turns per pole at load G, A-turn 
mean length of one field coil turn, in. 

used in armature conductor slot leakage permeance calculations 

total ampere turns, A-turn 

total ampere turns, A-turn 

main air gap ampere turns, A-turn 

main air gap ampere turns, A-turn 

maximum demagnetizing ampere turns per pole at rated load, A-turn 

direct-axis demagnetizing ampere turns per pole at load G, A-turn 

air-gap ampere turns per pole (no-load rated voltage), A-turn 

field current per coil, A 

NAMELIST name 

field current at rated load, A 

field winding resistance at temperature TF, ohm 

name of gas in alternator air gap 

used in distribution factor calculation 

pole ampere turns, A-turn 



+ FPLL 


pole ampere turns, A-turn 
FQ flux per pole (no-load rated voltage), kilolines 

+ FQLL flux per pole at no-load voltage QPERV, kilolines 

FRC rotor core ampere turns, A-turn 

+ FRCORL rotor core ampere turns, A-turn 

FS skew factor 

FS1 factor used in damper bar loss calculations 

FS2 factor used in damper bar loss calculations 

FSC short-circuit ampere turns for rated armature current, A-turn 

FTL tooth ampere turns, A-turn 

+ FTLL tooth ampere turns, A-turn 

*+ G volt-ampere alternator output at which load characteristics are calcu- 

lated, per unit or percent 

o 

GA main air-gap area, in. 

* GC minimum air gap (air gap at center of pole) (see fig. 12), in. 

GE effective main air-gap, in. 

GF constant used in load pole-face loss calculation 

GM factor used in load pole-face and damper loss calculations 

* GP maximum air gap (see fig. 12), in. 

GT ratio of stator slot opening width to main air-gap 

GXX flux per pole multiplying factor (see fig. 4) 

* H rectangular damper bar height, in. 

HC stator depth below slot, in. 

* HCOIL field coil height (see fig. 14), in. 

* HD damper bar slot opening height, in. 

HM armature conductor length (= 1/2 coil length), in. 

* HO stator slot dimension (see fig. 10), in. 

* HS stator slot dimension (see fig. 10), in. 

* HT stator slot dimension (see fig. 10), in. 

* HV number of stator cooling ducts 
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HW 

* HX 

* HY 
I 

IA 

IBN 

IC 

IDELR 

IDM 

IIQQ 

IPN 

* IPX 

* IQQ 

* IWF 

IZY 

IZZ 

J 

JA 

K 

KSAT 

KWIND 

LT 

* LTR 

* LTR1 

* LTS 

* M 
M 

MAGNET 

MM 


stator slot dimension (see fig. 10), in. 

stator slot dimension (see fig. 10), in. 

stator slot dimension (see fig. 10), in. 

used as index 

used as index 

number of damper bars 

number of parallel armature circuits per phase 

increment by which QPERV is increased, percent 

(integer) x (number of poles) x (number of phases) 

used to calculate distribution factor for fractional slot winding 

number of phases 

number of poles 

number of stator slots 

index to specify if windage loss is to be calculated (see table III under 
NAMELIST name AIRGAP) 

(number of poles) x (number of phases) 

type of stator slot (see fig. 10) 

used as index 

number of load calculations made before alternator saturates 
used as index 

saturation indicator (if KSAT = 10 alternator is saturated) 

index used in conjunction with subroutine OUTPUT 

dummy variable used in determination of stacking factors 

rotor pole body lamination thickness, in. 

rotor pole head lamination thickness, in. 

stator lamination thickness, in. 

molecular weight of gas in alternator air gap 

used as index 

subroutine name 

subscript of load point G such that G(MM) = 1.0 
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N 

NA 

NPF 

OUTPUT 

PI 

P2 

P3 

P4 

* PBA 

* PBH 

* PBL 

* PBW 
PC 
PDD 
PDQ 

* PE 
PEFF 
PF 
PF 

*+ PFC 

* PHH 

* PHL 

* PHW 
PI 
PN 
POL 


used as index 
used as index 

number of power factors for which load calculations are made 
subroutine name 

leakage permeance between inner pole head surfaces (see fig. 3), 
lines/A-turn 

leakage permeance between pole head end surfaces (see fig. 3), 
lines/A-turn 

leakage permeance between inner pole body surfaces (see fig. 3), 
lines/A-turn 

leakage permeance between pole body end surfaces (see fig. 3), 
lines/A-turn 

phase belt angle, deg 

pole body height (see fig. 13), in. 

pole body length (see fig. 13), in. 

pole body width (see fig. 13), in. 

specific armature conductor slot leakage permeance, lines/(A-turn)(in. ) 

specific permeance of damper bars in direct axis, lines/(A-turn)(in. ) 

specific permeance of damper bars in. quadrature axis, 
lines/(A-turn)(in. ) 

pole embrace 

effective pole leakage permeance = 2(P1 + P2) + P3 + P4, lines/A-turn 
design power factor 

power factor at which load characteristics are calculated 

power factors at which load characteristics are calculated 

pole head height (see fig. 13), in. 

pole head length (see fig. 13), in. 

pole head width (see fig. 13), in. 

rated line current, A 

number of phases 

dummy variable 
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+ PP pole -face losses under load, W 

PPL flux per pole in main air gap, stator core, and teeth, kilolines 

PPRL leakage flux per pole, kilolines 

+ PPRLL leakage flux per pole, kilolines 

+ PR field copper loss, W 

*+ PRESS array of gas pressures in alternator air gap corresponding to load 

2 

points G, lb/in. 

9 

PRSR gas pressure in alternator air gap, lb/in. 

2 

* PRSRAT gas pressure in alternator air gap at rated load VA, lb/in. 

+ PS armature copper loss, W 

* PT number of field coil turns per pole 

* PTCH stator winding pitch, per unit 

PX number of poles 

Q1 used in damper bar loss calculations 

Q2 used in damper bar loss calculations 

Q3 used in damper bar loss calculations 

Q4 used in damper bar loss calculations 

Q5 used in damper bar loss calculations 

QC used in Carter coefficient calculation 

QN slots per pole per phase 

+ QPERV voltage points on no-load saturation curve, percent 
QQ number of stator slots 

+ QVLL line-to-line voltage on no-load saturation curve, V rms 

+ QVLN line -to -neutral voltage on no-load saturation curve, V rms 

QZ used in damper bar loss calculations 

RATING NAMELIST name 

RB resistivity of stator conductor at temperature TTA, pohm-in. 

RC field coil weight (total of all coils), lb 

* RD field conductor diameter or width, in. 

* RE damper bar resistivity at 20° C, p.ohm-in. 
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RE1 

Reynolds number of gas in alternator air gap 

* 

RF 

type of armature winding (random or form wound) 


RG1 

armature winding resistance at temperature TST, ohm 

* 

RK 

pole body stacking factor 

* 

RK1 

pole head stacking factor 


RM 

damper bar resistivity at temperatures T3 and T33, pohm-in. 

*+ 

RMAT 

rotor material name 

* 

RO 

O 

density of gas in alternator air gap at rated load VA, lb/ft 


ROTOR 

NAMELIST name 

* 

RPM 

rotor speed, rpm 

* 

RR 

field coil resistivity at 20° C, pohm-in. 

+ 

RRA 

armature winding resistance, ohm 

+ 

RRB 

field coil resistance, ohm 

* 

RS 

armature conductor resistivity at 20° C, pohm in. 

* 

RT 

field conductor thickness, in. 


RY 

RRA/RB, Min." 1 


S 

2 

armature winding current density at rated load, A/in. 


SALENT 

main program name 


SALPHA 

sine of alpha 

* 

SB 

damper bar length, in. 

* 

SC 

number of armature conductors per stator slot 


SCR 

short circuit ratio 

* 

SD 

distance between centerlines of armature winding strands in depth (see 
fig. 11), in. 

* 

SEP 

clearance between armature winding end-turns (see fig. 11), in. 

* 

SF 

stator stacking factor 


SH 

armature conductor strand height (uninsulated), in. 


SI 

field self- inductance, H 


SIGMA 

2 

specific tangential force on rotor, lb/in. 

* 

SK 

stator slot skew at stator inside diameter, in. 
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SLOTS 

NAMELIST name 


SM 

tooth width at 1/3 distance from narrowest section, in. 

* + 

SMAT 

stator material name 

* 

SN 

strands per armature conductor in depth 

* 

SN1 

total strands per armature conductor 

+ 

SP 

total losses, W 


ss 

solid stator stack length, in. 

+ 

ST 

stator tooth loss, W 


STATET 

number of effective armature winding turns 


STATOR 

NAMELIST name 


STFK 

stacking factor (dummy variable) 


STRLEN 

length of straight part of armature winding end-turn (see fig. 18), in. 

* 

T1 

armature winding temperature at rated load, °C 

* 

Til 

armature winding temperature at no-load, °C 

* 

T2 

field winding temperature at rated load, °C 

* 

T22 

field winding temperature at no-load, °C 

* 

T3 

hot damper bar temperature, °C 

* 

T33 

cold damper bar temperature, °C 


TA 

armature time constant, sec 

* 

TB 

damper bar pitch, in. 


TC 

open-circuit time constant (field only), sec 

* 

TEMP 

temperature of gas in alternator air gap, °C 

* 

TF 

field coil temperature at which FK1 and TC are calculated, °C 


TG 

total air-gap flux at no-load, rated voltage, kilolines 


THETA 

angle used in pole weight calculations, rad 


TP 

pole pitch, in. 


TS 

stator slot pitch at stator inside diameter, in. 

* 

TST 

armature winding temperature at which RG1 and TA are calculated, 
°C 


TT 

stator slot pitch at 1/3 distance from narrowest section, in. 
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+ TTA 
+ TTB 
TTC 
U 

UA 

UU 

VI 

V2 

* VA 
VC 
VG 

* VIS 

* VMIN 

VR 
VS 
VT 
W 
W2 
W3 
W5 
+ WA 
WC 
WD 
WF1 
WFTTL 
WT 

WIND 

WINDGE 

WINDNG 


armature winding temperature, °C 
field winding temperature, °C 
transient time constant, sec 
used to calculate end-turn length 
= G(M), per unit 

used to calculate end-turn length 

used in damper bar loss calculations 

used in damper bar loss calculations 

rating' of alternator, kVA 

used in damper bar loss calculations 

used in damper bar loss calculations 

viscosity of gas in alternator air gap, lb/sec-ft 

minimum voltage at which no-load saturation curve is calculated, per 
unit 

rotor peripheral velocity at main air-gap, ft/min 

used in damper bar loss calculations 

used in damper bar loss calculations 

rating of alternator, W 

used in damper bar loss calculations 

used in damper bar loss calculations 

used in damper bar loss calculations 

alternator output power, kW 

armature winding weight (copper only), lb 

no-load damper loss at temperature T3, W 

windage loss in alternator air gap neglecting stator slots, W 

windage loss in alternator air gap including stator slots, W 

stator iron weight, lb 

NAMELIST name 

subroutine name 

NAMELIST name 
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WL 

WMIS 

WN 

WND 

WO 

WPOLE 

WQ 

WQL 

WRCORE 

WROTOR 

WT 

WTOTAL 

WU 

WW 

X 

XA 

XB 

XD 

XDD 

XDQ 

XF 

XI 

XL 

XN 

XQ 

XR 

XU 

YY 


core loss at flux density BK, W/lb 
miscellaneous load losses, W 
no-load pole-face loss, W 
windage loss, W 

damper bar slot opening width, in. 
weight of one rotor pole, lb 
no-load rated voltage core loss, W 
core loss, W 
weight of rotor core, lb 
rotor weight, lb 

no-load rated-voltage tooth loss, W 
total electromagnetic weight, lb 
no-load damper loss at temperature T33, W 
no-load damper loss (dummy variable), W 

flux density at which AT is found by interpolation (see fig. 17), 

2 

kilolines/in. 

synchronous reactance (direct), percent 
synchronous reactance (quadrature), percent 
armature reaction reactance (direct), percent 
damper leakage reactance (direct), percent 
damper leakage reactance (quadrature), percent 
field winding leakage reactance, percent 

angle between line current phasor and direct-axis (internal power fac- 
tor angle) (see fig. 4), rad 

stator winding leakage reactance, percent 

negative sequence reactance, percent 

armature reaction reactance (quadrature), percent 

reactance factor 

transient reactance (direct), percent 

subtransient reactance (direct), percent 



XX 

XXX 

XY 

Y 

* YY 
Z 

ZA 

ZB 

ZC 

ZD 

ZE 

ZF 

ZG 

ZY 

ZY 

* ZZ 


slope of magnetization curve at flux density X (see fig. 17) 
dummy variable used in variety of calculations 
subtransient reactance (quadrature), percent 
used in interpolation procedure of material curves 

slots spanned per armature coil (number of slots between coil sides +1) 
used in armature conductor slot leakage permeance calculations 
used in armature conductor slot leakage permeance calculations 
used in armature conductor slot leakage permeance calculations 
used in armature conductor slot leakage permeance calculations 
used in armature conductor slot leakage permeance calculations 
used in armature conductor slot leakage permeance calculations 
used in armature conductor slot leakage permeance calculations 
(length of field winding conductor )/(field conductor area), in. “* 

= 0.7HS, in. 

used as an index in armature and field conductor area calculations 
type of stator slot (see fig. 10) 
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TABLE I. - SUMMARY OF CALCULATIONS FOR PPL, FGL1, AND FGXL 



No-load rated- 
voltage magnetization 
characteristics 

Alternator load characteristic 
(rated voltage) 

No-load satura- 
tion data 

Flux per pole 

(PPL) 

PPL = FQ 

For power factor <0. 95 
PPL = FQ • (length of OP a ) 
For power factor >0. 95 
PPL = FQ • (length of 
OP) • 1. 10 

PPL = EDD ■ FQ 

Air-gap 
ampere-turns 
per pole 
(FGLl) 

FGLl = FH 

For power factor <0. 95 
FGLl = FH • (length of 
OP a ) 

For power factor >0. 95 
FGLl = FH • (length of 
OP) • 1. 10 

FGLl = EDD ■ FH 

Demagnetizing 
ampere -turns 
(direct axis) 
(FGXL) 

FGXL = 0 

FGXL = FGML • G • sin(XI) a 

FGXL = 0 


a OP and XI are obtained from vector diagram shown in fig. 4. 


TABLE II. - DEFINITION OF INPUT VARIABLES USED WITH MATERIAL DECKS 


Card 

FORTRAN symbol 

Format 

Description of data 

1 

SMAT 

6A6 

Stator material name 

2 - 5 

AI(J), 1<J<29 

8F10. 1 

AI(1) is the saturation flux density of the stator 
material; AI(2) to AI(29) are, alternately, 
values of flux density (kilolines/in. ) and 
magnetizing force (A-turn/in. ) for 14 points 
on the stator material magnetization curve; 
the 14 points are arranged in order of in- 
creasing flux density; note that AI(1) = AI(28) 

6 

RMAT 

6A6 

Rotor material name 

7 - 10 

AI(J), 3DsJ<59 

8F10. 1 

AI(3 1) is the saturation flux density of the rotor 

material; AI(32) to AI(59) are, alternately, 

o 

values of flux density (kilolines/in. ) and 
magnetizing force (A-turn/in. ) for 14 points 
on the material magnetization curve; the 
14 points are arranged in order of increasing 
flux density; note that AI(31) = AI(58) 












TABLE III. - DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

name 

FORTRAN 

symbol 

Definition 

Classifi- 

cation 

(a) 

Remarks 

RATING 

VA 

Kilovolt-ampere rating of alternator, kVA 

M 




EE 

Line-to-line design voltage, rms V 

C 


Either 

one must be read in, or both may be 


EP 

Line-to-neutral design, voltage, rms V 

C 

1 

read in 


F 

Frequency, Hz 

C' 





RPM 

Shaft rotational speed, rpm 

C 


Any two must be read in, or all three may 


IPX 

Number of poles 

Q, 

J 

be read in 


G 

Volt -ampere load at which load character- 

0 

G is a 

subscripted variable (array size is 5); 



istics are calculated, percent or per unit 


if not read in, program assumes values, 





0, 0.75, 1.0, 1.25, and 1. 50; any one or 





all (except 0) may be changed by reading 





in different values; program automatically 





arranges values in increasing order; any 





number >9. 0 is assumed to be in percent, 





s9. 

0 in per unit. If, at some value of G, 





the alternator saturates, then G will be 





decreased in steps of 0. 1 until saturation 





no longer occurs 


PFC 

Power factors at which load character- 

M 

PFC Is a subscripted variable (array size 



istics are calculated 


is 5); one or more values must be read 





in; 

the first value will be assumed to be j 





the design power factor 


VMIN 

Lowest voltage for which a point on the open- 

o 

If not read in, VMIN = 0. 70 will be assumed 



circuit saturation curve will be calculated, 






per unit 




STATOR 

DI 

Stator inside diameter, in. 





DU 

Stator lamination outside diameter, in. 





CL 

Length of stator stack, in. 

1 




HV 

Number of cooling ducts 

mm 

If there are no cooling ducts, these need 


BV 

Width of cooling duct, in. 

■31 

not be read in 


SF 

Stacking factor (stator) 

c 

■ 

Either one or both may be read in; if 


LTS 

Stator lamination thickness, in. 

c 

■ 

neither is read in, program assumes 





■ 

that stator is not laminated (SF = 1. 0) 


WL 

Core loss at flux density BK, W/lb 





BK 

Flux density at which core loss WL is 






given, kiloltnes/in. ^ 




SLOTS 

zz 

Stator slot type 

M 

See fig. 10 


BO 

Slot dimension (use with slot types 2, 3, 

C 





and 4), in. 






B3 

Slot dimension (use with slot type 3), in. 






BS 

Slot dimension (use with slot types 1, 2, 







4, and 5), in. 






HO 

Slot dimension (use with slot types 2, 3, 







and 4), In. 






HX 

Slot dimension (use with slot types 1, 2, 







3, and 5), in. 






HY 

Slot dimension (use with slot types 1, 2, 







3, and 5), in. 






HS 

Slot dimension (use with slot types 1, 2, 3, 

M 





4, and 5), in. 





HT 

Slot dimension (use with slot types 2 and 3), 

C 




IQQ 

Number of stator slots 

M 




a M, mandatory; C, conditional; O, optional. 
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TABLE III. - Continued. DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

FORTRAN 

Definition 

Classifi- 

Remarks 

name 

symbol 


cation 





(a) 


WINDNG 

RF 

Type of armature coil 

M 

RF =1.0 for form wound coil; RF = 0 for 





random wound coil 


SC 

Number of conductors per stator slot 

M 

In fig. 11, SC = 4 


YY 

Slots spanned per coil (number of slots 
between coil sides plus one) 

I 

Either one must be read in, or both may 


PTCH 

Armature winding pitch, per unit 

cj 

be read in 



Number of parallel circuits per phase 

M 







DW 

Strand diameter or width, in. 

M 

See fig. 11 


SN 

Strands per armature conductor in depth 

C 

Read for rectangular wire only (In fig. 11, 



(radial direction) 


SN = 2) 


SN1 

Total strands per armature conductor 

M 

In fig. 11, SN = 4 


DW1 

Uninsulated strand thickness (radial di- 

C 

Read for rectangular wire only; see fig. 11 



rection), in. 




CE 

Armature coil extension, in. 

M 

See fig. 11 


SD 

Distance between centerline of strands in 

M 

See fig. 11 



depth, in. 




PBA 

Phase belt angle, deg 

0 

If not read in, program assumes PBA - 60° 


SK 

Stator slot skew at stator inside diameter, 
in. 

Hated-load armature winding temperature, 

o 

If not read in, program assumes SK = 0 


T1 

M 

Used for loss and efficiency calculations 


RS 

Armature conductor resistivity at 20° C, 

O 

If not read in, program assumes copper re- 



(/iohm)(in. ) 


sistivity (0. 694) 


ALPHAS 

Armature conductor temperature coef- 

o 

If not read in, program assumes copper 



ficient of resistivity at 20° C, (°C)~* 


temperature coefficient (0. 00393) 


Til 

No-load armature winding temperature, 

M 

Used for loss and efficiency calculations 


TST 

Armature winding temperature, °C 

0 

Program calculates and prints out armature 





resistance at this temperature; if not read 
in, program assumes TST = 25° C 


EL 

Armature coil end turn length, in. 

0 

Read in if exact value is known; if not, pro- 





gram will calculate approximate value 


BX 

Thickness of armature coil at end-turn, 

o 

See fig. 11; if not read in, program as- 



in. 


sumes BX = BS - 0. 015 


SEP 

Clearance between armature winding end 

o 

See fig. 11; if not read in, program as- 



turns, in. 


sumes SEP = 0. 01 

AIRGAP 

GC 

Minimum air gap (at center of pole), 
in. 

Maximum air gap, in. 

M 

See fig. 12 


GP 

c 

Need not be read in if air gap is constant 





(i. e. , if GP = GC); see fig. 12 


IWF 

Index to specify if windage loss is to be 

o 

To have program calculate approximate 



calculated 


windage loss in main air gap, set IWF = 1 
and add data cards for air-gap gas name 
and NAMELIST name WIND at end of de- 
sign deck 


a M, mandatory; C, conditional; O, optional. 
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TABLE m. - Continued. DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

name 

FORTRAN 

symbol 

Definition 

Classifi- 

cation 

(a) 

Remarks 

ROTOR 

RK 

Stacking factor of pole body 

C 1 


One or both may be read in; if neither is 


LTR 

Lamination thickness of pole body, in. 

c i 


read in program assumes that pole body 




J 


is not laminated (RK =1.0) 


RK1 

Pole head stacking factor 

> 

c 

1 

One or both may be read in; if neither is 


LTR1 

Pole head lamination thickness, in. 

c 


read in, program assumes solid pole 






head (RK1 = 1. 0) 


PBL 

Pole body length (axial direction), in. 

ci 


If PBL = PHL, only one (either one) need be 


PHL 

Pole head length (axial direction), in. 

c 


read in; see fig. 13 


PE 

Pole embrace 

cl 


One must be read in; both may be read in; 


PHW 

Pole head width, in. 

C J 


see fig. 13 


PHH 

Pole head height, in. 

M 

See fig. 13 


PBH 

Pole body height, in. 

M 

See fig. 13 


PBW 

Pole body width, in. 

M 

See fig. 13 


DISH 

Inside rotor diameter (for hollow rotor), 

C 

Read in only for rotor with center hole (see 



in. 


fig. 14) 


DSH 

Outside diameter of rotor core, in. 

O 

Program calculates value if not read in 





(see fig. 14) 


WROTOR 

Rotor weight, lb 

O 

If not read in, program will calculate ap- 





proximate rotor weight 


D1 

Pole face loss factor 

0 

If not read in, D1 is calculated from value of 





LTR1 using the following: D1 = 1. 17 for 





LTRl 5 0. 045; D1 = 1. 75 for 0. 045 < 





LTR1 £ 0. 094; D1 = 3. 5 for 0. 094 < 





LTRl £ 0. 17; D1 = 7. 0 for LTRl > 0. 17; 





if LTRl is not read in, program calcu- 





lates value of LTRl based on RK1 


Cl 

Ratio of fundamental maximum to actual 

o'' 





maximum value of field form (field 






form is air-gap flux density distribu- 



Constants defined for conventional salient 



tion due to field only) 



pole alternator (ref. 13) if not an input, 


CP 

Ratio of average to maximum value of 

O 

m 

values are calculated from formulas 



field form 


m 

given in refs. 4 and 13 


CM 

Demagnetizing factor (direct axis) 

0 




CQ 

Cross magnetizing factor (quadrature axis) 

■ 

i 


FIELD 

HCOIL 

Field coil height, in. 

0 

If not read in, program assumes 





HCOIL = PBH, see fig. 14 


PT 

Number of field coil turns per pole 

M 



RD 

Field conductor diameter or width, in. 

M 



RT 

Field conductor thickness, in. 

c 

Do not read in for round conductors 


BCOIL 

Field coil width, in. 

M 

See fig. 14 


T2 

Rated-load field temperature, °C 

M 


Used in loss and efficiency calculations 


T22 

No-load field temperature, °C 

M 

i 



RR 

Field-coil resistivity at 20° C, 

0 

If not read in, 0. 694 is assumed 



(pohm)(in. ) 




ALPHAR 

Temperature coefficient of resistivity at 

o 

If not read in, 0. 00393 is assumed 



20° C, “C' 1 




TF 

Field-coil temperature, °C 

0 

Program calculates and prints out field- 





coil resistance and open-circuit time 





constant at this temperature; if not read 





in, program assumes TF = 25° C 


FE 

Mean length of one turn, in. 

0 

Program calculates approximate value if 





not read in 


a M, mandatory; C, conditional; O, optional. 
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TABLE III. - Concluded. DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

name 

FORTRAN 

symbol 

Definition 

Classifi- 

cation 

(a) 

Remarks 

DAMPER 

BN 

Number of damper bars per pole 

M 

If BN = 0, none of following variables for 






DAMPER need be read in 


WO 

Damper bar slot opening width, in. 

C 




HD 

Damper bar slot opening height, in. 

C 




DD 

Damper bar diameter, in. 

C 


For round damper bars only 


H 

Rectangular damper bar height, in. 

C 

} 

For rectangular damper bars only 


B 

Rectangular damper bar width, in. 

c 

J 



SB 

Damper bar length, in. 

c 


If not read in, SB = PHL is assumed 


TB 

Damper bar pitch, in. 

c 


If not read in, approximate value is calcu- 






lated 


T33 

Cold damper bar temperature, °C 

c 


If not read in, 20° C will be assumed 


T3 

Hot damper bar temperature, °C 

c 




RE 

Damper bar resistivity at 20° C, 

c 


If not read in, program assumes copper re- 



(pohm)(in. ) 



sistivity (0. 694) 


ALPHAE 

Temperature coefficient of resistivity 

c 


If not read in, program assumes copper 



at 20° C, °C' 1 



temperature coefficient (0. 00393) 

WIND* 3 

VIS 

Viscosity of gas in alternator air gap, 

M 




lb/sec-ft 





M 

Molecular weight of gas 





TEMP 

Temperature of gas, °C 





PRSRAT 

Pressure of gas at alternator rating 






VA, lb/in. 2 





PRESS 

Array of gas pressures corresponding 

o 

If not read in, program scales PRSRAT di- 



to alternator load points G, lb/in. 2 



rectly with load 


RO 

Density of gas at alternator rating VA, 

0 

If not read in, program calculates value 



lb/ft 3 





a M, mandatory; C, conditional, O, optional. 

^This card is needed only if IWF of NAMELIST name AIRGAP is unity. It must be preceded by a card giving the gas name 


(see p. 13). 


Fjrj 















Figure 1. - Simplified flow chart of salient, wound pole alternator computer program. 
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Figure 8. - Material deck for SAE 4340 steel. 


Points whose coordinates 
are used in preparation 
of material deck (see 
table for exact values) 

■ Magnetization curve 


Magnetic flux Magnetizing 


density 

force 

a 0 

a 0. 01 

6.45 

12.0 

12.90 

24.0 

19.35 

29.0 

58.0 

42.0 

64.5 

44.0 

77.4 

50.0 

83.8 

56.0 

90.3 

70.0 

96.8 

92.0 

103.0 

126.0 

109.5 

184 

116.0 

280 

122.5 

454 


Maximum \ 
magnetic 
flux density 


“Not shown on graph 


4 6 8 10 20 40 60 80 100 200 

Magnetizing force, ampere-turns/in. 

Figure 9. - Average magnetization curve for SAE 4340 steel. 









HY ~i 

i 

r~i 

k 

HX 

□ 

i 

HS 



□ 

i 


BS 



BO* 

JB_S = Bi 


Types land 5: Open 
slot, constant slot 
width. Type 5 slot is 
same as type 1, but 
it contains only one 
coil side. 



Type 2: Partly closed slot, 
constant slot width. 



HO 


T-7 


1 P" 80 


f 

HS 

I 


BS-*- 



Type 3: Partly closed slot, Type 4: Round slot, 

constant tooth width. 


Figure 10. - Stator slot dimensions. (Starred variables are not 
input; they are shown for reference only. ) 



Figure 11. - Definition of variables used with NAMELIST nameWINDNG. 
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Damper bar 
slot opening 7 



Damper bar length 

5.20 

Damper bar pitch 

0. 213 

Number of damper bars per pole 

9 

Damper bar material 

Copper 

Hot damper bar temperature, °C 

150 

Cold damper bar temperature, °C 

20 


(a) Damper bar design parameters. 


BN 9 

WO 0.030 

HD 0,070 

DD 0. 100 

H Not read in when damper bars are round 

B Not read in when damper bars are round 

SB 5.20 

TB 0.213 

T33 Need not be read in ; program automatically assumes 20° C 
T3 150 

RE Need not be read in; program automatically assumes value for copper 

ALPHAE Need not be read in ; program automatically assumes value for copper 



SDAH 

p 

CR BH=9, UD=0.03> HD=0. 0?» 10=0.10. SB=5.e(l, TB=0.813. T3=130 S 

- — 


C-^.2 

a 





s 

FORTRAN STATEMENT 
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(b) Conversion of design parameters to FORTRAN symbols and data card $DAMPER. 
Figure 15. - Preparation of data card for NAMELIST name DAMPER. 


$DATA 

VANADIUM pERMENOUR 


154. 


12.9 

1.92 

38.7 


2.62 


77.4 


3.23 


90.3 


3.53 


1C3. 

4.35 

109.7 


5.25 


116. 


& . 66 


122. 

5 

8.68 


129. 

12.5 

135.5 


20.2 


142. 


44.4 


145. 

3 

131. 


148.3 

363. 

154. 


2320. 








SAE 4340 

STEEL IRGCM 

TEMP) 











122.5 


O. 

O.Ol 

6.45 


12. 


12.90 


24. 


19. 

35 

29. 


58. 

42. 

64.5 


44. 


77.4 


50. 


83. 

8 

56. 


90.3 

70. 

96.8 


92. 


103. 


126. 


109 

. 5 

184. 


116. 

28C. 

122.5 


454. 








tRATING 

VA= 

200, EP=170» F=500 » 

I PX=l 2 

» G 

=0, .5, 1. 

C, 1 

.5,2.0, 

PFC 

= .75, 1 

.0 



$ S T A TOR 

D I = 

12.5, DU= 

14.5, CL=5, 

,20, SF 

=0. 

96, L T b = 

.006, WL = 1 1 

.8, 

BK = 7 7 

. 4 

$ 


tSLCTS 

11 = 2 

, BO=0 .07 

, BS=0.142, 

, H0=0 • 

03, 

HX=0 . 345 , 

HY=0 . 10, 

HS 

=0.455 

, HT=0 . 

05! 

IQQ 

= 180 

S 












SWINONG 

RF = 

l, SC=2 , 

C = 3, SN= l , 

SN 1 = 1 » 

DW 

=0.112, 

DW 1 

=3.150, 

CE = 

<■'•25, 

SD 

-o 

o 

ii 

0 , 

Tl = 

20C, 

T1 1= 100, 

PTCH=3.8 

S 











tAIRGAP GOC.07, IWF = 1 S 

!RCTGR 01=1. 17, PHL=5.20, PE=C.6667 t DISH=8.2, DSH=9.9, PHH=^.4&, PBH=0.75, 
PBW= 1 . 30 » RK1=0.9C, LTR1=.0^6 $ 

IF I ELD BCCI L =0 • 42 , PT=44, RO = 0.l9, RT=0.025, T2=200, T22=100 l 
fOAMPER BN = 9, W0=0.03» HC=0.37, DD=0.iO, SB=5.20, TB=0.213, T3=L50 t 
AIR 

$ W I NO M=2 9 . 9 i V l S= 1 . 2 1 E-05 * D RSRAT=15.0, TEMP=150 * 


Figure 16. - Listing of data deck for sample alternator design. 



Typical straight line segment 
between input data points on 
material magnetization curves-. 



I 

/ 

/ 

/ 

/ 

AI(NA)— 

/ 

\ 

x-\ 


ANNA-2H 

7 

Slope 

XX 



/ 

/ 

/ 



AI(NA-l) 


AT 

AKNA+l) 


Log of magnetizing force, ampere-turns/in. 


Figure 17. - Variables used in interpolation 
between points on magnetization curves. 



length calculations. 
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